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INTRODUCTION 

Cancer  patients  with  no  clinical  evidence  of  disease  after  the  initial  treatment  frequently  relapse 
with  distant  metastasis  years  later.  Prior  to  diagnosis  and  treatment,  primary  tumors  may  release 
large  numbers  of  cancer  cells  into  the  circulation.  Although  a  majority  of  the  dispersed  cells 
perish  in  the  bloodstream  or  soon  after  infiltrating  distant  organs,  a  minority  may  survive  as 
latent  seeds  in  host  tissues  (7).  As  a  result,  cancer  patients  who  are  clinically  considered  disease- 
free  carry  thousands  of  disseminated  tumor  cells  (DTCs)  in  the  bone  marrow  and  other  organs(2, 
3).  DTCs  are  frequently  negative  for  proliferation  markers,  suggesting  that  they  survive  in  a 
quiescent  state  that  is  resistant  to  cytotoxic  adjuvant  treatments^,  5).  Little  is  known  about  the 
nature  of  dormant  DTCs  and  the  mechanisms  that  allow  them  to  remain  quiescent,  evade 
immunity,  and  retain  tumor-initiating  capacity.  Latent  metastasis  therefore  remains  a  largely 
unknown  biological  process  as  well  as  a  major  concern  in  the  clinic,  and  very  particularly  in 
breast  cancer. 

Our  ability  to  address  these  questions  and  to  gain  an  urgently  needed  understanding  of  the 
molecular  basis  for  latent  metastasis  have  been  limited  by  a  scarcity  of  preclinical  models  that 
recapitulate  key  features  of  this  metastatic  stage,  including  a  propensity  of  DTCs  to  enter 
proliferative  quiescence,  to  survive  long-term  as  quiescent  DTCs,  and  to  have  a  finite  probability 
of  initiating  macrometastatic  growth  ( 6 ,  7).  The  goal  of  this  innovative  award  is  to  address  this 
important  unmet  need  by  dissecting  and  targeting  the  biology  of  latent  metastasis. 

KEYWORDS 

Breast  cancer;  latent  metastasis;  disseminated  cancer  cells;  circulating  cancer  cells;  cancer  stem 
cells;  metastatic  niche;  SOX  genes;  WNT  pathway;  NK  lymphocytes. 

OVERALL  PROJECT  SUMMARY 

Progress  under  Task  1:  To  develop  experimental  models  of  latent  metastasis. 

Progress  Item  1.1.  Latency  competent  cells  isolated  form  early-stage  breast  and  lung 
cancers.  We  isolated  cancer  cells  that  are  competent  to  seed  relevant  organs  with  latent 
metastasis  (latency  competent  cancer  cells,  LCC  cells).  As  sources,  we  used  cell  lines  derived 
from  early-stage  human  breast  and  lung  adenocarcinomas.  The  H2087  line  was  established  from 
a  stage  I  lung  aden o c arc  i  n om a( 8) ,  and  the  HCC1954  line  from  a  stage  IIA  HER2+  breast  cancer 
(9).  Nearly  half  of  early-diagnosed  (stage  I  and  II)  lung  adenocarcinoma  patients  develop  distant 
relapse  despite  surgical  resection  of  the  primary  tumor,  implying  the  presence  of  latent 
disseminated  disease  (10).  The  HER2+  breast  cancer  patient  population  is  experiencing  a  marked 
increase  in  the  incidence  of  brain  metastasis  after  anti-HER2  targeted  therapies  that  effectively 
suppress  extracranial  relapse  and  extend  survival  (77).  Thus,  both  cancer  types  are  important 
sources  of  latent  metastasis  in  the  clinic. 

H2087  and  HCC1954  cells  transduced  with  a  GFP-luciferase  reporter  and  antibiotic  resistance 
vectors  were  intracardially  injected  into  the  arterial  circulation  of  FoxnI'm  mice.  The  mutation  in 
Foxnl  renders  the  mice  athymic,  severely  blunting  the  maturation  of  effector  T  cells  but 
preserving  innate  immunity  components  including  macrophages  and  NK  cells  (72).  As 
monitored  by  bioluminescence  imaging  (BLI),  the  luciferase  signal  became  undetectable  7  days 
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after  inoculation,  indicating  extensive  clearance  of  the  cancer  cells  from  the  circulation  or  after 
organ  infiltration.  Most  mice  remained  signal-free  and  healthy  for  3  months.  Organs  from  these 
mice  were  dissociated  into  single-cell  suspensions  in  culture  to  recover  antibiotic-resistant  cancer 
cells.  We  were  able  to  recover  these  cells  from  lung  (H2087-LCC1)  and  kidney  (H2087-LCC2) 
of  mice  injected  with  H2087,  and  from  brain  of  mice  injected  with  HCC1954  (HCC1954-LCC1). 
The  tumorigenic  activity  of  LCC  cells  after  orthotopic  implantation  was  similar  to  that  of  the 
parental  populations.  Thus,  LCC  cells  retained  tumor  initiating  potential  after  months  of  latency 
in  mouse  tissues. 

Of  20  athymic  mice  injected  with  H2087-LCC  cells,  one  developed  overt  metastasis,  one 
developed  a  spinal  metastasis  after  4  months  and  did  not  progress,  two  developed  incipient 
lesions  after  7  months,  and  16  remained  free  of  detectable  BLI  signal  for  over  8  months  (Figure 
1A,C).  Only  1  out  of  8  mice  injected  with  HCC1954-LCC1  developed  overt  metastasis,  and  the 
rest  remained  metastasis-free  for  over  4  months  (Figure  1B-C).  In  comparison,  aggressive 
metastatic  lines  that  we  previously  isolated  from  late-stage  human  breast  carcinoma  (MDA-MB- 
231)  and  lung  adenocarcinoma  (H2030)  formed  extensive  metastases  within  3  weeks  in  all  mice 
(13,  14). 

We  were  able  to  histologically  detect  and  recover  cancer  cells  from  20  of  26  (77%)  LCC-injected 
mice  that  remained  metastasis-free,  but  only  from  2/14  (14%)  of  mice  injected  with  the  parental 
lines.  H2087-LCC1  and  -LCC2  targeted  both  the  lungs  and  the  kidneys,  suggesting  that  the  latent 
phenotype  of  these  cells  was  not  strictly  organ  specific  (Figure  1D-G).  Cells  isolated  from  rare 
macrometastatic  tumors  arising  in  LCC-injected  mice  (H2087-LCC-M  and  HCC1954-LCC-M 
lines)  showed  a  latent  phenotype  upon  reinjection,  with  no  increase  in  macrometastatic  activity, 
arguing  that  LCCs  stochastically  form  overt  metastases. 

Progress  Item  1.2.  LCC  cell  localization  in  infiltrated  organs.  We  chose  H2087-LCC1, 
H2087-LCC2,  HCC1954-LCC1,  and  the  HCC1954-LCC2  line  derived  from  LCC1  injected 
mice,  for  further  analysis.  Two  weeks  after  inoculation,  approximately  20%  of  disseminated 
GFP+  HCC1954  LCC  cells  in  the  brain  were  found  in  small  clusters  (>10  cells)  and  the  rest  as 
single  cells.  The  proportion  of  cell  clusters  remained  stable  or  progressively  decreased  over  time. 
After  3  months,  more  than  90%  of  the  DTCs  were  single  cells  and  the  rest  were  small  clusters  of 
approximately  20  cells,  as  determined  by  anti-GFP  immunostaining  (Figure  1D-E),  H&E 
staining  (Figure  IF),  and  human  Vimentin  immunostaining  (Figure  1G).  In  the  kidney,  LCC  cells 
were  detected  between  the  renal  tubules,  adjacent  to  capillaries,  or  within  glomeruli  (Figure  1F- 
G).  A  few  LCC  cells  were  observed  within  the  adrenal  gland  cortex.  In  the  lung,  LCC  cells  were 
invariably  found  within  the  alveolar  walls  (Figure  1G).  In  the  brain  and  other  organs,  LCC  cells 
were  closely  associated  with  the  vasculature  (Figure  ID).  Aggressive  brain  metastatic  cells 
rapidly  spread  over  the  abluminal  surface  of  microcapillaries,  and  this  spreading  is  required  for 
colony  outgrowth  (15).  In  contrast,  LCC  cells  infiltrating  the  brain  parenchyma  spread  along 
capillaries  only  transiently  (days  1-3  after  inoculation),  and  subsequently  adopted  a  rounded 
morphology  (Figure  1H-I). 

Progress  Item  1.3.  LCC  cell  entry  into  quiescence.  To  monitor  proliferation  in  the  first  days 
after  inoculation,  we  labeled  LCC  cells  with  5-ethynyl-2'-deoxyuridine  (EdU)  prior  to 
intravenous  injection  into  athymic  mice  (Figure  2A-B).  Proliferation  would  dilute  the  amount  of 
EdU  retained  in  these  cells.  After  14  days,  approximately  60%  of  LCC  cells  that  reached  the 
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lungs  still  retained  EdU  versus  15-25%  in  the  parental  populations  (Figure  2C-D),  indicating  that 
LCC  cells  entered  quiescence  more  readily.  Three  months  after  inoculation,  approximately  90% 
of  the  H2087-LCC1  cells  in  lung  and  HCC1954-LCC1  cells  in  brain  were  negative  for  the 
proliferation  marker  Ki-67.  Ki-67+  LCC  cells  detected  were  largely  confined  to  cell  clusters 
(Figure  2E-F).  To  determine  the  effect  of  a  mitogen-poor  environment  we  cultured  LCC  cells  in 
mitogen- low  media  (MLM,  2%  serum)  or  regular  mitogen-rich  media  (MRM,  10%  serum). 
Under  MLM  conditions,  LCC  cells  underwent  a  rapid  decrease  in  proliferation,  as  determined  by 
dye  retention  (Figures  2G),  whereas  the  parental  line  showed  little  (H2087)  or  no  decrease 
(HCC1954)  in  proliferation.  LCC  cells  accumulated  with  a  G0/G1  DNA  content  (Figure  2H)  and 
no  change  in  apoptosis  marker  (cleaved  caspase-3)  levels. 

Progress  Item  1.4.  Stem  cell  nature  of  LCC  cells.  Gene  Set  Enrichment  Analysis  of 
differentially  expressed  genes  in  HCC1954-LCC  cells  compared  to  parental  HCC1954  identified 
a  mammary  stem  cell  signature  (MaSC)  as  a  top-scoring  gene  set,  under  both  MLM  and  MRM 
conditions.  Principal  component  analysis  revealed  strong  clustering  of  HCC1954-LCC  cells  with 
human  and  mouse  mammary  stem  and  progenitor  cells,  away  from  mature  luminal  and  stromal 
compartments  (Figure  3 A)  (16).  Parental  HCC1954  clustered  with  the  mature  luminal 
compartment.  Moreover,  HCC1954-LCC  cells  were  enriched  for  the  surface  marker  profile 
CD44Hl/CD24Lo,  which  is  typical  of  human  breast  cancer  stem  cells  (Figure  3B)  (17).  Based  on 
gene  expression  profiles  associated  with  molecularly  distinct  cell  types  of  the  mouse  lung 
epithelium  (18),  H2087-LCC  cells  clustered  with  the  stem-like  alveolar  type  I  (ATI)  and 
bipotent  progenitor  (BP)  cells.  Parental  H2087  clustered  with  alveolar  type  II  (AT2)  and  Clara 
cell  lineages  under  MRM  and  MLM  conditions,  respectively,  (Figure  3C).  These  results 
suggested  that  our  LCC  isolation  protocol  enriched  for  cancer  cells  with  stem/progenitor  cell 
markers. 

In  sum,  LCC  cells  recapitulate  key  features  of  the  latent  metastatic  state,  including  a 
propensity  to  enter  proliferative  quiescence,  an  ability  to  survive  as  latent  entities  in  relevant 
organs  for  months,  and  to  retain  tumorigenic  capacity  and  metastasis-initiating  potential. 

Progress  under  Task  2:  Identify  genes  that  are  specifically  active  in  latent  metastatic  cells. 

Progress  Item  2.1.  Sox2  and  Sox9  association  with  the  LCC  phenotype.  Progenitor  cell 
identity  is  determined  by  lineage-specific  transcription  factors  (19).  Notably,  two  master 
regulators  of  stem  and  progenitor  cell  identity,  Sox2  (20)  and  Sox9  (21),  ranked  high  among 
transcription  factors  whose  expression  was  prominently  associated  with  the  LCC  phenotype. 
Sox2  was  predominant  in  H2087-LCC  cells  whereas  HCC1954-LCC  cells  showed  high  Sox9 
expression  and  low  Sox2  expression  (Figures  3D-E).  Sox2  and  Sox9  regulate  stem  and 
progenitor  cells  in  adult  tissues  (19)  and  Sox2  mediates  pluripotency  reprogramming  in 
differentiated  cells  (22).  Both  genes  are  implicated  in  cancer.  Sox2  is  genetically  amplified  in 
squamous  carcinomas  and  small  cell  lung  cancers  (23,  24).  Sox9  is  upregulated  in  brain  tumors 
and  basal  cell  carcinomas  (25),  and  confers  stem  cell  phenotype  to  breast  cancer  cells  (21). 
Analysis  of  TCGA  datasets  confirmed  a  high  incidence  of  Sox2  and  Sox9  genetic  alterations 
(amplifications)  in  different  tumor  types,  including  breast  and  lung  cancers. 

Progress  Item  2.2.  Metastatic  seeding  by  LCC  cells  requires  Sox  factors.  We  generated 
H2087  and  HCC1954-LCC  cell  lines  expressing  shRNAs  targeting  Sox 2  and  Sox9  respectively. 
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Sox2  knockdown  in  H2087-LCC  and  Sox9  knockdown  in  HCC1954-LCC  cells  did  not  alter  the 
growth  rate  of  these  cells  in  monolayer  culture.  We  inoculated  H2087-LCC  cells  into  athymic 
mice  by  tail  vein  injection,  and  monitored  mice  by  whole  body  BLI  imaging  weekly  for  60  days. 
No  BLI  signal  was  detected  during  this  time  period  or  on  sacrificing  the  mice  and  imaging  lungs 
ex  vivo.  Further  analysis  of  the  lungs  showed  a  marked  decrease  in  the  metastatic  seeding 
capacity  of  Sox2  depleted  H2087-LCC  cells  compared  to  cells  expressing  a  control  vector 
(Figure  3F).  Similarly,  Sox9  depleted  HCC1954-LCC1  cells  showed  a  marked  decrease  in  the 
ability  to  seed  the  brain  (Figure  3G).  The  ability  of  tumor-initiating  cells  to  form  colonies 
(oncospheres)  in  suspension  provides  a  test  of  their  capacity  to  initiate  tumor  growth  under 
restrictive  conditions  (26).  Parental  H2087  and  HCC1954  cells  generated  the  same  number  of 
oncospheres  as  their  LCC  counterparts.  Interestingly,  Sox2  expression  was  enriched  in  H2087- 
LCC  oncospheres,  but  not  oncospheres  formed  by  parental  H2087  (Figure  3H).  Sox2  knockdown 
inhibited  oncosphere  formation  in  H2087-LCC  but  not  in  parental  (Figure  31).  Thus,  LCC  cells 
are  distinct  subpopulations  of  stem-like  cancer  cells  characterized  by  high  expression  of  Sox2  or 
Sox9  and  a  dependence  on  these  transcription  factors  for  growth  under  restrictive  conditions  in 
vitro  and  for  seeding  of  latent  metastasis  in  mice. 

Progress  under  Task  3:  Identification  of  key  stromal  components  in  latent  metastasis. 

Progress  Item  3.1.  Metastatic  latency  balanced  by  NK  cell  immune  surveillance.  Although 
LCC  cells  were  superior  at  seeding  various  organs,  they  still  suffered  massive  attrition  in 
athymic  mice.  The  elimination  of  disseminated  cancer  cells  may  be  due  to  metabolic  and 
mechanical  stresses  (6,  7),  but  also  to  immune  surveillance  (27-29).  Foxnlnu  athymic  mice  retain 
functional  NK  cells  and  other  components  of  the  innate  immune  system,  raising  the  possibility 
that  immune  surveillance  restricts  the  expansion  of  LCC  cells  in  these  mice.  To  test  this 
possibility,  we  inoculated  LCC  cells  into  NOD.Cg -Prkdcscid  Il2rgtmIWjl/ SzJ  mice  (NOD/SCID 
Gamma,  NSG  mice),  which  are  defective  for  both  adaptive  and  innate  immune  responses  (30). 
Strikingly,  H2087-LCC  and  HCC1954-LCC  cells  gave  rise  to  overt  metastasis  with  high 
penetrance  (Figure  4A-B)  and  in  multiple  organs  in  NSG  mice  (Figure  4C-D).  These  organs 
included  the  liver,  a  site  normally  rich  in  NK  cells  (31).  To  test  the  role  of  NK  cells  as  candidate 
inhibitors  of  LCC  expansion,  we  depleted  athymic  mice  of  NK  cells  by  administration  of 
polyclonal  anti-asialo-GMl  antibody  or  the  anti-NKl.l  monoclonal  antibody  PK136  (32,  33). 
Both  NK  cell  depletion  strategies  resulted  in  permissive  outgrowth  of  LCC  cells,  as  manifested 
by  a  marked  increase  in  whole  body  BLI  signal,  and  in  the  ability  to  grow  in  the  brain  as  shown 
by  ex  vivo  brain  imaging  and  quantification  of  average  brain  lesion  size  (Figures  4E-F).  Similar 
results  were  observed  upon  NK  cell  depletion  in  mice  injected  with  H2087-LCC2,  which 
resulted  in  bone  metastasis  (Figures  4G). 

Progress  Item  3.2.  NK  cell-enforced  latency  in  syngeneic  models.  We  also  investigated  the 
role  of  NK  cells  as  suppressors  of  DTC  expansion  in  syngeneic  immunocompetent  mouse 
models.  One  model  was  based  on  lung  adenocarcinoma  cell  lines  derived  from  tumor-bearing 
KrasG12D/p53Del  mice  (34).  The  482T1  (T-Met)  cell  line  readily  metastasizes  to  the  liver  when 
injected  into  the  spleen  of  syngeneic  B6129SF1  mice,  whereas  the  368T1  (T-nonMet)  line  does 
not  form  overt  metastases  in  these  assays  (34).  We  conducted  these  experiments  in  control 
B6129SF1  mice  or  mice  that  were  treated  with  anti-asialo-GMl  antibody  to  deplete  NK  cells. 
NK  cell  depletion  resulted  in  a  100-fold  increase  in  the  liver  metastatic  activity  of  T-nonMet 
cells,  reaching  a  level  that  was  comparable  to  that  of  T-Met  cells  (Figure  4H-I).  We  obtained 


6 


similar  results  in  4T07  cells,  a  non-metastatic  cell  line  derived  from  a  spontaneous  BALB/c 
mouse  mammary  tumor  (35). 

Collectively,  these  results  suggested  that  LCC  cells  can  proliferate  after  infiltrating  distant 
organs,  but  NK  cell  immune  surveillance  prevents  the  accumulation  of  progeny  and  only  spares 
LCC  cells  that  enter  quiescence. 

Progress  under  Task  4:  Establish  the  functional  relevance  of  LCC  genes  and  pathways. 

Progress  Item  4.1.  Downregulation  of  NK  cell  activators  in  quiescent  LCC  cells.  A  balance 
of  activating  and  inhibitory  signals  regulates  the  ability  of  NK  cells  to  target  cancer  cells  for 
elimination  (36,  37).  Oncogenic  transformation  leads  to  loss  of  NK  cell  inhibitory  receptors  and 
upregulation  of  NK  cell  activating  ligands  in  tumor  cells,  which  in  turn  promotes  NK  cell- 
mediated  clearance  of  cancer  cells.  During  tumor  progression,  cancer  cell  evolve  to  escape  NK 
cell  mediated  recognition  (38-41).  As  LCC  cells  survived  NK  cell-mediated  immune  surveillance 
in  athymic  mice,  we  queried  LCC  gene  expression  profiles  for  potential  mechanisms  of  immune 
evasion. 

We  examined  gene  expression  signatures  associated  with  immune  recognition  and  clearance  by 
macrophages,  T,  B,  or  NK  cells.  Genes  signatures  associated  with  NK  cell-mediated 
cytotoxicity,  in  particular,  were  downregulated  in  LCC  cells  compared  to  parental  H2087  and 
HCC1954  populations.  Moreover,  quiescence-inducing  (MLM)  conditions  induced  additional 
changes  in  the  expression  of  NK  cell  activators  in  LCC  cells.  The  genes  that  were  further 
downregulated  in  quiescent  LCC  cells  encode  key  NK  cell  activators  (Figure  5A-B).  These  genes 
include  UL  16-binding  proteins  (ULBP;  also  known  as  retinoic  acid  early  transcript,  RAET) 
ULBP2/RAET 1 H,  ULBP3/RAET1N,  and  ULBP5/RAET 1 G,  are  frequently  expressed  by 
malignant  cells  and  stimulate  anti-tumor  responses  mediated  by  the  NK  activating  receptor 
NKG2D/CD314  (42),  and  PVR/CD  155  is  a  ligand  for  the  NK  cell  receptor  CD226/DNAM-1 
(DNAX  accessory  molecule- 1)  that  activates  cancer  cell  killing  (43).  The  pro-apoptotic  cytokine 
receptors  FAS  and  TRAILR,  which  are  critical  for  NK  cell-mediated  target  killing  (44,  45),  were 
also  downregulated  in  LCC  cells  under  quiescence  conditions  (Figure  5A-B).  To  determine  if 
LCC  cells  entering  quiescence  are  intrinsically  resistant  to  the  cytotoxic  action  of  NK  cells,  we 
incubated  parental  and  LCC  cells  in  low-mitogen  media  and  then  added  freshly  isolated,  IL-2 
activated  mouse  spleen  NK  cells  to  the  cultures  (Figure  5C).  Compared  to  parental  populations, 
LCC  cells  showed  resistance  to  cytolysis  when  incubated  with  NK  cells  (Figure  5D).  Thus,  LCC 
cells  that  enter  quiescence  undergo  a  striking  downregulation  of  a  set  of  NK  cell  activators,  and 
acquire  resistance  to  NK  mediated  cell  killing. 

Progress  Item  4.2.  Attenuated  WNT  signaling  associated  with  LCC  quiescence.  To  identify 
what  primes  LCC  cells  to  enter  a  quiescent  state,  we  applied  signaling  pathway  classifier 
analysis  to  the  transcriptomic  data  sets  (46).  Under  quiescence  conditions,  all  H2087-LCC  cell 
lines  showed  a  marked  reduction  in  WNT,  MYC  and  NF-kB  signaling  activity,  and  an  increase 
in  TGF-P  signaling  activity  (Figure  6A).  WNT  signaling  response  was  also  attenuated  in 
HCC1954  LCC  cells  under  these  conditions.  MYC  is  a  downstream  target  of  WNT  signaling 
(47),  and  TGF-P  is  a  known  mediator  of  cytostasis  (48).  A  high  ratio  of  phospho-p38  to 
phospho-ERK  kinases,  previously  described  in  dormant  cancer  cell  models  (49),  was  present  in 
HCC 1 954-LCC  but  not  in  H2087-LCC  cells. 
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WNT  is  a  potent  mitogen  for  stem  and  progenitor  cells,  and  is  implicated  in  metastatic  outgrowth 
in  lung  adenocarcinoma  (14)  and  breast  cancer  stem  cells  (50,  51).  H2087-LCC  cells 
demonstrated  a  clear  resistance  to  pathway  activation  by  WNT3A  addition,  as  determined  by 
Axin2  expression  (Figure  6B)  (52),  and  TCF  transcriptional  reporter  activation  (53).  Therefore, 
we  investigated  the  possibility  that  LCC  cells  produce  a  WNT  inhibitor.  We  observed  that  the 
WNT  signaling  inhibitor  dickkopf-related  protein  1  (DKK1)  was  highly  expressed  in  LCC  cell 
lines  compared  to  parental  counterparts,  both  at  the  mRNA  and  protein  levels  (Figure  6C). 
DKK1  expression  by  LCC  cells  in  vivo  was  confirmed  by  immunofluorescence  staining  (Figure 
6D). 

Progress  under  Task  5:  Establish  the  clinical  relevance  of  LMBC  genes  and  pathways 

Progress  Item  5.1.  We  have  started  this  task  with  work  on  a  set  of  breast  cancer  patient  derived 
xenografts  (PDXs),  which  recapitulate  the  histopathology,  heterogeneity  and  metastatic 
properties  of  the  primary  tumor  (54,  55).  Surgical  orthotopic  implantation  of  hormone  receptor¬ 
negative  breast  cancers  PDX  models  HCI-001  and  HCI-002,  and  HER2+  breast  cancer  PDX 
model  HCI-008,  gave  rise  to  tumors  with  no  visible  metastasis  (Figure  6E)  (54).  Orthotopic  PDX 
tumors  expressed  human  Vimentin,  which  we  used  as  a  surrogate  marker  to  identify 
disseminated  cancer  cells  in  the  lungs,  brain  and  kidneys  of  mice  bearing  these  tumors.  14%  - 
51%  of  disseminated  cancer  cells  were  positive  for  both  human  vimentin  and  DKK1  expression 
(Figure  6E).  These  data  demonstrated  DKK1  expression  in  disseminated  tumor  cells  from 
orthotopically  implanted  PDX  models. 

Progress  under  Task  6:  Mechanism  of  action  of  LCC  genes  and  pathways. 

Progress  Item  6.1.  Direct  mechanistic  link  between  SOX2  and  DKK1  expression.  DKK1  is  a 
direct  target  of  Sox2 -mediated  transcriptional  activation  in  mesenchymal  stem  cells,  where  Sox2 
binds  to  the  ~  75  base  pairs  upstream  of  the  DKK1  transcriptional  start  site  (56,  57).  Indeed, 
H2087-LCC1  cells  showed  Sox2  binding  to  this  promoter  region  (Figure  6F).  Stable  shRNA 
knockdown  of  Sox2  in  LCC  cells  decreased  the  expression  of  DKK1  (Figure  6G).  DKK1 
knockdown  re-sensitized  LCC  cells  to  WNT  pathway  activation,  as  demonstrated  by  increased 
Axin2  expression  upon  WNT3A  addition  (Figure  6H).  Dye  retention  assays  demonstrated  that 
DKK1  knockdown  stimulated  the  proliferation  of  LCC  cells  under  low-mitogen  conditions 
(Figure  61). 

Progress  Item  6.2.  Autocrine  DKK1  enforces  a  quiescent,  immune  evasive  state  in  latent 
metastasis.  The  preceding  results  suggested  that  autocrine  DKK1  primes  LCC  cells  to  enter 
quiescence,  which  leads  to  downregulation  of  NK  cell  activators.  One  prediction  of  this  model  is 
that  autocrine  DKK1  is  required  to  protect  LCC  cells  from  elimination  in  athymic  mice.  In  line 
with  this  prediction,  shRNA-mediated  DKK1  depletion  in  LCC  cells  significantly  decreased  the 
accumulation  of  LCC  cells  in  the  lungs  of  inoculated  athymic  mice  (Figure  7A).  Moreover,  Ki- 
67+  stromal  cells  surrounded  the  cancer  cells  in  bronchiolar  regions  under  these  conditions. 
Many  of  these  stromal  cells  were  positive  for  the  leukocyte  common  antigen  CD45,  suggesting 
that  DKK1 -depleted  LCC  cells  trigger  lymphocyte  infiltration  and  for  cancer  cell  clearance. 
DKK1  knockdown  in  H2087-LCC  cells  increased  the  expression  of  NK  cell  activating  ligands 
ULBP1,  2,  4  and  5,  and  several  death  signal  receptors  (Figure  7B).  Moreover,  depletion  of  NK 
cells  in  athymic  mice  rescued  the  formation  of  macrometastases  by  DKK1 -knockdown  H2087- 
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LCC2  cells  (Figure  7C),  suggesting  that  LCC  cells  liberated  from  DKK1  growth  inhibition  can 
grow  as  macrometastases  in  a  growth  permissive  microenvironment. 

In  contrast  to  the  deleterious  effect  of  DKK1  knockdown  on  the  viability  of  LCC  cells  in  athymic 
mice,  DKK1  knockdown  increased  the  metastatic  growth  of  LCC  cells  in  NSG  mice  (Figure  7D- 
E).  Conversely,  overexpressing  DKK1  in  H2087-LCC  cells  inhibited  their  growth  in  NSG  mice. 
As  a  corollary,  depletion  of  Sox2  in  H2087-LCC  or  Sox9  in  HCC1954-LCC  cells  inhibited  the 
ability  of  these  cells  to  form  metastasis  in  NSG  mice  (Figure  7F)  as  it  did  in  athymic  mice.  These 
data  confirmed  the  essential  role  of  Sox  transcription  factors  for  metastatic  growth  of  LCC  cells, 
and  the  ability  of  NK  cells  to  restrain  LCC  outgrowth  but  spare  quiescent  LCC  cells  (Figure  7G). 

Progress  under  Task  7:  Define  the  feasibility  of  targeting  residual  disease  genes  and 
pathways  in  pre-clinical  trials 

Two  lines  of  work  that  we  conducted  in  parallel  with  the  above  studies,  and  partly  supported 
with  this  Award,  have  begun  to  address  the  feasibility  of  targeting  different  aspects  of  residual 
metastatic  disease  based  on  our  basic  research  findings. 

Progress  Item  7.1.  Pharmacologic  inhibition  of  gap  junction  activity  against  brain 
metastasis.  Brain  metastasis  represents  a  substantial  source  of  morbidity  and  mortality  in  various 
cancers,  and  is  characterized  by  high  resistance  to  chemotherapy.  The  role  of  the  most  abundant 
cell  type  in  the  brain,  the  astrocyte,  in  brain  metastasis  has  remained  ambiguous.  In  ongoing 
studies  we  have  demonstrated  that  breast  cancer  and  lung  cancer  cells  express  protocadherin  7 
(PCDH7)  to  favor  the  assembly  of  carcinoma-astrocyte  gap  junctions  composed  of  connexin  43 
(Cx43).  Once  engaged  with  the  astrocyte  gap-junctional  network,  brain  metastatic  cancer  cells 
employ  these  channels  to  pass  cytosolic  dsDNA  to  the  astrocytes,  activating  the  cGAS-STING 
pathway  and  production  of  inflammatory  cytokines  IFNa  and  TNFa.  As  paracrine  signals,  these 
factors  activate  the  STAT1  and  NF-kB  pathways  in  brain  metastatic  cells,  which  support  tumor 
growth  and  chemoresistance. 

The  evidence  that  genetic  inhibition  of  gap  junction  components  decreased  brain  metastatic 
outgrowth  provided  a  rationale  for  testing  pharmacologic  suppressors  of  gap  junction  activity 
against  brain  metastasis.  To  this  end,  we  selected  two  orally  bioavailable  compounds  for  pre- 
clinical  trials.  In  addition  to  anti-inflammatory  activity,  meclofenamate  inhibits  Cx43  gap 
junction  gating  (55),  inhibits  epileptogenesis  in  animal  models  (59),  passes  the  BBB  after 
systemic  administration  (55),  is  well  tolerated  systemically  (60)  and  is  currently  an  FDA- 
approved  NSAID.  Tonabersat  is  an  benzopyran  derivative  that  binds  to  a  unique  stereoselective 
binding  site  in  astrocytes  (61,  62),  inhibits  gap-junction-mediated  pathophysiological  processes 
including  cortical  spreading  depression  (63)  and  trigeminal  ganglion  neuronal-satellite  cell 
signaling  in  animal  models  (64),  and  was  systemically  well-tolerated  and  safe  in  patients  with 
migraine  (65). 

Both  Tonabersat  and  meclofenamate  inhibited  dye  transfer  from  astrocytes  to  cancer  cells  as 
measured  by  flow  cytometry  (Fig.  8A),  and  the  release  of  IFNa  and  TNFa  in  co-cultures  of  these 
cells  (Fig.  8B),  recapitulating  the  phenotype  seen  in  knockdown  of  Cx43  or  PCDH7.  Mice  were 
treated  with  either  vehicle  or  with  these  compounds  from  day  1  following  arterial  inoculation  of 
MDA231-BrM2  cells  or  H2030-BrM3  cells  in  immunodeficient  mice,  or  KRas/p53-393Nl  cells 
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in  immunocompetent  mice  (Fig.  8C).  Both  drugs  prevented  the  emergence  of  brain  metastases, 
consistent  with  our  evidence  that  gap  junction  activity  is  relevant  for  metastatic  outgrowth. 
However,  this  treatment  did  not  restrict  growth  of  MDA231-BrM2  cells  as  lung  metastatic 
lesions  or  as  orthotopic  tumors. 

Progress  Item  7.2.  Gap  junction  directed  therapy.  To  test  the  effect  of  Cx43  or  PCDH7 
depletion  in  established  metastases,  we  transduced  MDA231-BrM2  cells  with  Tet-inducible 
shRNA  expression  vectors  (Fig.  8E).  A  red  fluorescence  protein  (RFP)  under  the  control  of  the 
same  promoter  provided  a  marker  of  hairpin  expression  in  vivo.  Cells  transduced  with  inducible 
Cx43  or  PCDH7  shRNA  vectors  showed  doxycycline-dependent  depletion  of  Cx43  or  PCDH7, 
respectively.  These  cells  were  injected  intracardially  and  allowed  to  form  brain  metastases  for  14 
days.  At  this  stage,  brain  lesions  are  apparent  by  BLI  in  all  mice;  the  aggressive  lesions  engulf 
the  microvasculature  (Fig.  8D)  and  will  result  in  death  of  the  animals  in  2-3  weeks  (75). 
Doxycycline  administration  starting  on  day  14  resulted  in  reduced  brain  metastatic  burden  three 
weeks  later,  compared  to  controls  (Fig.  8F,G). 

Brain  metastases  are  distinguished  by  pronounced  resistance  to  chemotherapy  (66,  67). 
Carboplatin  crosses  the  BBB  (68),  with  modest  improvement  in  overall  survival  in  patients  with 
brain  metastases  from  breast  (69)  or  lung  cancer  (70).  Carboplatin  alone  (50  mg/kg/5  days) 
starting  on  day  14  inhibited  brain  metastasis  to  a  similar  extent  as  depletion  of  Cx43  or  PCDH7 
(Fig.  8F,G);  combination  carboplatin  and  doxycycline  reduced  the  metastatic  burden  further  (Fig. 
8F,G).  Therefore,  we  assessed  the  effectiveness  of  combination  gap  junction  modulatory  therapy 
with  chemotherapy  (Fig.  8H).  Treatment  with  carboplatin  alone  minimally  inhibited  brain 
metastasis  growth  (Fig.  81).  Either  Tonabersat  (10  mg/kg)  or  meclofenamate  (20  mg/kg)  as 
single  agents  (Fig.  81)  significantly  inhibited  progression  of  metastatic  lesions  at  the  35-day  end 
point.  The  combination  of  carboplatin  with  either  Tonabersat  or  meclofenamate  profoundly 
inhibited  brain  metastasis  (Fig.  81). 

Progress  Item  7.3.  Pilot  Clinical  Trial  NCT  02429570.  With  these  results  we  introduce  a 
therapeutic  strategy  based  on  existing  modulators  of  gap  junctions,  meclofenamate  and 
Tonabersat,  that  breaks  this  paracrine  loop,  and  provide  proof-of-principle  for  its  applicability  to 
treat  established  brain  metastasis.  Based  on  these  results,  the  co-first  author  of  the  submitted 
manuscript,  Dr.  Adrienne  Boire,  MD,  PhD,  and  a  neurononcologist,  has  initiated  a  Pilot  Clinical 
Trial  (NCT  02429570)  of  Meclofenamate  in  subjects  with  recurrent  or  progressive  brain 
metastasis  from  solid  tumor  primary  https://clinicaltrials.gov/ct2/show/NCT02429570 

Progress  Item  7.4.  Therapy-induced  tumor  secretomes  and  resistance  and  tumor 
progression.  Drug  resistance  invariably  limits  the  clinical  efficacy  of  targeted  therapy  with 
kinase  inhibitors  against  cancer.  In  recently  published  work  (77)  we  show  that  targeted  therapy 
with  BRAF,  ALK,  or  EGFR  kinase  inhibitors  induces  a  complex  network  of  secreted  signals  in 
drug-stressed  melanoma  and  lung  adenocarcinoma  cells.  This  therapy-induced  secretome  (TIS) 
stimulates  the  outgrowth,  dissemination,  and  metastasis  of  drug-resistant  cancer  cell  clones  and 
supports  the  survival  of  drug-sensitive  cancer  cells,  contributing  to  incomplete  tumour 
regression.  The  vemurafenib  reactive  secretome  in  melanoma  is  driven  by  down-regulation  of 
the  transcription  factor  FRA1.  In  situ  transcriptome  analysis  of  drug-resistant  melanoma  cells 
responding  to  the  regressing  tumour  microenvironment  revealed  hyperactivation  of  multiple 
signalling  pathways,  most  prominently  the  AKT  pathway.  Dual  inhibition  of  RAF  and 
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PI3K/AKT/mTOR  pathways  blunted  the  outgrowth  of  the  drug-resistant  cell  population  in  BRAF 
mutant  melanoma  tumours,  suggesting  this  combination  therapy  as  a  strategy  against  tumour 
relapse.  Thus,  therapeutic  inhibition  of  oncogenic  drivers  induces  vast  secretome  changes  in 
drug-sensitive  cancer  cells,  paradoxically  establishing  a  tumour  microenvironment  that  supports 
the  expansion  of  drug-resistant  clones,  but  is  susceptible  to  combination  therapy.  (Refer  to 
Appendix  -  Publication). 

KEY  RESEARCH  ACCOMPLISHMENTS 

•  Selection  for  phenotypically  representative  models  of  latent  metastasis  from  early-stage 
lung  and  breast  cancers,  Latency  Competent  Cancer  (LCC)  cells 

•  LCC  cells  shown  to  have  stem-like  properties,  are  prone  to  enter  quiescence,  and  highly 
express  Sox2  or  Sox9  stem/progenitor  identity  factors.  Disseminated  cancer  cells  with 
these  characteristics  can  persist  long-term  as  latent  seeds  of  metastasis. 

•  Demonstration  that  LCC  cells  can  initiate  metastasis  after  infiltrating  distant  organs 

•  Demonstration  that  NK  cell-mediated  immune  surveillance  prevents  the  accumulation  of 
LCC  progeny,  but  spare  LCC  cells  that  enter  quiescence. 

•  Discovery  that  LCC  cells  express  DKK1  as  an  autocrine  promoter  of  cell  quiescence. 

•  Finding  that  LCC  quiescence  leads  to  downregulation  of  NK  cell  activators  and  evasion 
of  NK  cell-mediated  killing. 

•  Proof  of  principle  that  gap  junction  modulators  tonabersat  and  meclofenamate  can  inhibit 
brain  metastasis 

•  Initiation  of  Clinical  Trial  NCT  02429570,  Meclofenamate  in  subjects  with  recurrent  or 
progressive  brain  metastasis  from  solid  tumor  primary 

CONCLUSION 

The  overall  goals  of  this  project  are  to  biologically  dissect,  molecularly  deconstruct,  and 
conceptually  understand  latent  metastasis  of  breast  cancer,  in  order  to  target  its  viability  for  the 
prevention  of  metastasis.  Our  progress  towards  this  goal  is  on  schedule.  We  have  accomplished 
the  tasks  that  were  scheduled  for  the  period  of  this  Progress  Report.  We  have  encountered  new 
and  unexpected  findings.  There  have  been  no  major  unforeseen  problems  in  the  development  of 
this  work. 

CHANGES  IN  PERSONNEL 

•  Saloni  Agrawal  has  been  added  to  the  project  as  a  Research  Technician 
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metastasis  through  the  cytosolic  dsDNA  response  pathway.  Submitted  for 
publication 

3.  Invited  Articles: 

•  Obenauf,  A.  and  Massague,  J.  Surviving  at  a  distance:  organ  specific  metastasis. 
Trends  in  Cancer,  in  press  NIHMS714178 

4.  Abstracts:  Nothing  to  report 

b.  Presentations  made  during  the  last  year. 

b.l.  International: 

•  Charles  Rodolphe  Brupbacher  Symposium,  Breakthroughs  in  Cancer  Research  and 
Therapy ,  Zurich  Switzerland 

•  Gordon  Research  Conference  on  Stem  Cells  &  Cancer,  Ventura,  CA 

•  ISPEN  Cancer  Series  Meeting,  Tumor  Heterogeneity  and  Tumor  Microenvironment, 
Mysore,  India 

•  EMBO/EMBL  Symposium,  Frontiers  in  Stem  Cells  &  Cancer,  Heidelberg,  Germany 

•  AACR  Conference  on  Advances  in  Brain  Cancer  Research,  Washington,  DC 
b.2.  National: 

•  Origins  of  Metastatic  Traits,  Danny  Thomas  Lecture,  St.  Jude  Children’s  Hospital, 
Memphis,  TN 

•  Deconstructing  Metastasis,  Distinguished  Guest  Lecture,  Baylor  College  of  Medicine, 
Houston,  TX 

•  Deconstructing  CNS  Metastasis,  Neurology  Grand  Rounds,  Albert  Einstein  College  of 
Medicine,  NY 

•  Scientific  Perspectives  on  Cancer,  Surgery  Grand  Rounds,  MSKCC,  New  York,  NY 

•  Stock  Lecture,  Medicine  Grand  Rounds,  MSKCC,  New  York,  NY 
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•  Metastasis  Initiating  Cells,  Cell  Press  Lab  Links  Intersections  of  Stem  Cell  Biology  and 
Cancer,  NYU,  NY 

INVENTIONS,  PATENTS  AND  LICENSES 

•  Massague,  J,  Boire,  A.,  and  Chen,  Q.,  Methods  For  Treating  Brain  Metastasis.  United 
States  Provisional  Patent  Application  No.  62/052,966  Filed  September  19,  2014. 

REPORTABLE  OUTCOMES 

•  Identification  of  DKK1  and  NK  activators  as  potential  therapeutic  targets  to  eliminate 
disseminated  cancer  cells. 

•  Proof  of  principle  that  gap  junction  modulators  tonabersat  and  meclofenamate  can  inhibit 
brain  metastasis,  and  initiation  of  a  clinical  trial. 

OTHER  ACHIEVEMENTS 

Joan  Massague: 

2014  National  Prize  for  Research  in  Biology,  Spain 

2015  Charles  Rodolphe  Brupbacher  Prize  for  Cancer  Research 

2015  Keynote  address,  Stem  Cells  and  Cancer  Gordon  Conference 

2015  Member,  External  Advisory  Board,  Salk  Institute  Cancer  Center 
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APPENDICES 

Figure  1.  Localization  and  dormancy  of  disseminated  LCC  cells. 

(A)  BLI  tracking  of  mice  injected  with  H2087-LCC1,  H2087-LCC2,  or  the  aggressive  metastatic 
lung  cancer  line  H2030-BrM  for  comparison.  Each  line  represents  an  individual  mouse  tracked 
over  the  indicated  time  period.  Mice  with  BLI  signal  above  the  dotted  red  line  (whole  body 
photon  flux  average  radiance  >106  photons/sec/cm2/sr)  were  euthanized. 

(B)  BLI  tracking  of  mice  injected  with  HCC1954-LCC1  over  5  months,  including  the  aggressive 
metastatic  breast  cancer  line  MDA231-BrM  for  comparison. 

(C)  Representative  images  of  BLI  signal  in  mice  injected  with  the  indicated  LCC  cell  lines.  The 
number  of  mice  represented  by  each  image  is  indicated  at  the  bottom. 

(D)  Representative  immunofluorescence  staining  of  disseminated  LCC  cells  (green)  next  to 
capillaries  (red)  in  the  parenchyma  of  BLI-negative  organ  sites. 

(E)  H2087-LCC2  cells  (green)  infrequently  detected  as  cell  clusters  in  the  lung  (top  panel)  or 
kidney  (bottom  panel)  remain  associated  with  the  vasculature  (red). 

(F)  Spatial  localization  of  single  H2087-LCC1  cells  (arrows)  in  BLI-negative  organs  3  months 
post-injection  by  H&E  staining.  LCC  cells  were  stochastically  distributed  in  the  renal  tubule 
(top-left)  or  glomeruli  (top-right)  of  the  kidney,  within  alveolar  walls  of  the  lung  (bottom-left),  or 
capillaries  in  the  brain  (bottom-right). 

(G)  Immunohistochemical  staining  of  disseminated  H2087-LCC1  cells  from  panel  IF  with 
human  vimentin  (brown)  to  confirm  the  presence  of  human  cancer  cells  in  these  BLI-negative 
organs. 

(H)  Time-course  analysis  of  HCC1954-LCC1  cell  morphology  (green)  as  they  extravasate  and 
colonize  the  brain.  LCC  cells  were  intimately  associated  with  brain  capillaries  (red),  initially 
spread  along  capillaries  (left  inset)  followed  by  quick  adoption  of  a  rounded  phenotype  (right- 
inset). 

(I)  Quantification  of  HCC1954-LCC1  cell  morphology  in  the  brain  3  days  and  14  days  post¬ 
injection.  Data  are  mean  percentage  of  coopting  cells  per  brain  ±  S.E.M.  N  =  3  mice  per  group, 
scoring  representative  serial  sections  of  the  entire  brain  for  each  mouse.  P  <  0.001  (***), 
Student’s  t-test.  Scale  bars,  10pm  (D),  100pm  (E),  20pm  (F  and  G),  and  50pm  (H). 


•  Figure  2.  LCC  cells  adopt  a  slow-cycling  state  in  vitro  and  in  vivo. 

(A)  Experimental  design  for  EdU  pulse-chase  experiment.  LCC  cells  were  labeled  with  EdU  for 
24hr  in  vitro,  and  injected  into  mice  by  tail-vein  injection.  Lungs  were  harvested  after  2  weeks, 
and  representative  serial  sections  for  each  lung  were  quantitated  for  EdU+  LCC  cells. 

(B)  LCC  cells  (green)  efficiently  labeled  in  vitro  with  EdU  (red). 

(C)  Detection  of  double-positive  EdU+/GFP+  LCC  cells  versus  parental  counterparts  in  lungs 
harvested  2  weeks  post-injection.  Data  are  mean  percentage  of  EdU+  cells  per  lung  ±  S.E.M.  N  = 
3  mice  per  group,  scoring  representative  serial  sections  of  the  entire  lung  for  each  mouse. 
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(D)  Representative  immunofluorescence  images  of  double-positive  GFP+/EdU+  LCC  cells 
(green,  red)  quantified  in  Figure  2C.  P  <  0.05  (*),  P  <0.01  (**),  Mann-Whitney  Test. 

(E)  Representative  immunofluorescence  images  of  proliferating  Ki-67+  (red)  H2087-LCC1  (in 
lung)  and  HCC1954-LCC1  cells  (in  brain)  3  months  post-injection. 

(F)  Quantification  of  LCC  cells  from  Figure  2E.  Data  are  percentage  of  total  cells  that  are  Ki67+ 
per  organ  ±  S.E.M.  N  =  5  per  group,  scoring  representative  serial  sections  of  the  entire  organ  of 
each  mouse.  P  <  0.001  (***),  Student’s  t-test. 

(G)  Retention  of  eFluor670  dye  by  H2087  and  HCC1954  derivatives  after  6  days  of  in  vitro 
culture  in  MRM  or  MLM  conditions. 

(H)  Cell  cycle  analysis  of  H2087  and  HCC1954  derivatives  by  BrdU/APC  after  3  days  of  in  vitro 
culture  under  MRM  or  MLM  conditions.  Scale  bars,  10pm  (B  and  D),  and  15pm  (E). 


•  Figure  3.  LCC  cells  are  enriched  for  stem  cell-like  characteristics. 

(A)  Principal  Component  Analysis  (PC A)  of  HCC1954  derivatives  and  normal  breast  cell 
populations  of  human  or  mouse  origin  ( 16)  plotted  in  the  same  space.  HCC1954-LCC  cells 
resemble  mammary  stem  cell  (MaSC)  gene  expression  profiles. 

(B)  Flow  cytometry  analysis  of  cell  surface  CD24  and  CD44  expression  in  HCC1954  parental 
and  LCC1  cells.  HCC1954-LCC1  cells  show  marked  enrichment  in  the  CD44Hl/CD24Lo 
compartment,  indicative  of  breast  cancer  stem  cells. 

(C)  PCA  of  H2087  derivatives  and  mouse  bronchiolar  and  alveolar  cell  lineages  (Treutlein, 
2014).  H2087-LCC  cells  cluster  with  the  stem-like  alveolar  type  1  (AT)  and  bipotent  progenitor 
(BP)  cells. 

(D)  Sox2  and  Sox9  mRNA  expression  levels  in  H2087  and  HCC1954  derivatives. 

(E)  Sox2  and  Sox9  protein  expression  levels  in  H2087  and  HCC1954  derivatives. 

(F)  Sox2  depletion  attenuates  survival  of  H2087-LCC  cells  in  athymic  mice.  Data  are  total 
number  of  LCC  cells  scored  in  the  lung  3  months  post-injection,  with  whisker  plots  representing 
minimum  and  maximum  values.  N  =  5-6  mice  per  group,  scoring  representative  serial  sections  of 
the  entire  lung  for  each  mouse.  P  <0.001  (***),  Mann-Whitney  Test. 

(G)  Sox9  depletion  attenuates  survival  of  HCC1954-LCC1  cells  in  athymic  mice.  Data  are  total 
number  of  LCC  cells  scored  in  the  brain  2  months  post-injection,  with  whisker  plots  representing 
minimum  and  maximum  values.  N=  5  mice  per  group,  scoring  representative  serial  sections  of 
the  entire  brain  for  each  mouse.  P  <0.05  (*),  Mann-Whitney  Test. 

(H)  Sox2  expression  (red;  white  arrows)  in  oncospheres  generated  from  H2087  Parental  (top 
panel)  or  LCC1  (lower  panel)  cells. 

(I)  Oncosphere-forming  capacity  is  upon  Sox2  depletion  in  H2087-LCC  derivatives.  Side  panel 
depicts  representative  control  or  Sox2-depleted  H2087-LCC  cell  oncospheres.  P  <  0.01  (**), 
Student’s  t-test.  Scale  bars,  25pm  (H)  and  20pm  (I). 
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•  Figure  4.  LCC  cells  are  in  equilibrium  with  the  innate  immune  system. 

(A-B)  BLI  tracking  of  NSG  mice  injected  with  HCC1954-LCC  (A)  or  H2087-LCC  (B)  cells, 
with  aggressive  metastatic  lines  MDA231-LM2  or  H2030-BrM  shown  for  comparison.  Each  line 
represents  one  mouse  tracked  over  the  indicated  time  period. 

(C)  Representative  images  of  BLI  signal  in  NSG  mice  injected  with  the  indicated  LCC  cell  lines 
over  the  indicated  time  period. 

(D)  Metastatic  outgrowth  of  HCC1954-LCC1  (green)  in  the  liver,  brain,  kidney  and  lungs  of 
NSG  mice  2  months  post-injection.  Vasculature  visualized  in  red  with  isolectin-B4. 

(E)  Depletion  of  NK  cells  by  either  anti-NKl.l  (PK136)  or  anti-asialo-GMl  antibody  regimen 
allows  outgrowth  of  HCC1954-LCC1  cells  injected  in  athymic  nude  mice.  Outgrowth  quantified 
by  whole  body  photon  flux  at  2  months  post- injection,  with  whisker  plots  representing  minimum 
and  maximum  values.  N  =  5-7  mice  per  group.  P  <0.001  (***),  Mann-Whitney  Test. 

(F)  Representative  ex  vivo  brain  BLI  images  of  mice  from  Figure  4E.  Frozen  brain  tissue 
sections  from  these  mice  show  massive  outgrowth  of  HCC1954-LCC1  cells  (green)  upon  NK 
depletion. 

(G)  NK  cell  depletion  regimens  also  allow  for  outgrowth  of  H2087-LCC2  cells  in  athymic  nude 
mice.  Outgrowth  quantified  by  whole  body  photon  flux  at  2  months  post- injection,  with  whisker 
plots  representing  minimum  and  maximum  values.  N  =  5-7  mice  per  group.  P  <0.05  (*),  P  <0.01 
(**),  Mann-Whitney  Test. 

(H)  NK  cell  depletion  by  anti-asialo  GM1  allows  liver  outgrowth  of  the  non-metastatic  cell  line 
T-nonMet  intra-splenetically  injected  into  syngeneic  B6129SF1  mice.  The  metastatic  line  T-Met 
is  shown  for  comparison.  Outgrowth  quantified  by  whole  body  photon  flux  at  1  month  post¬ 
injection,  with  whisker  plots  representing  minimum  and  maximum  values.  N  =  5-7  mice  per 
group,  P  <0.01  (**),  Mann-Whitney  Test. 

(I)  Representative  ex  vivo  gross  (top  panel)  and  BLI  (bottom  panel)  liver  images  from  mice  in 
Figure  4H.  Scale  bars,  100pm  (D  and  F). 


•  Figure  5.  LCC  cells  evade  NK  cell  mediated  immune  surveillance. 

(A-B)  A  diverse  set  of  genes  important  for  NK  cell  recognition  and  cytotoxicity  are 
downregulated  in  both  H2087-LCC1  (A)  and  HCC1954-LCC1  (B)  when  grown  in  MLM 
conditions. 

(C)  Schematic  of  in  vitro  NK  cell  cytotoxicity  experiment.  Labeled  cancer  cells  grown  in  MLM 
conditions  are  incubated  with  activated  primary  NK  cells  at  a  1:5  target: effector  ratio  for  3  hours, 
followed  by  quantification  of  remaining  intact  cancer  cells. 

(D)  H2087  and  HCC1954-LCC  cells  are  more  resistant  to  NK-cell  mediated  cytotoxicity 
compared  to  parental  counterparts.  Data  are  mean  percentage  of  cytolysis  ±  S.E.M.  for  three 
replicates.  P  <0.01  (**),  P  <0.001  (***),  Student’s  t-test. 
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•  Figure  6.  Attenuated  WNT  signaling  enforces  quiescence  in  LCC  cells. 

(A)  Signaling  pathway  response  signature  analysis  of  H2087  derivatives  under  MLM  conditions 
shows  marked  reduction  in  WNT  and  MYC  pathway  activity  (upper  panel).  MYC  mRNA  z  score 
of  these  derivatives  was  scored  from  RNA-seq  data  (lower  panel). 

(B)  H2087-LCC  cells  are  less  responsive  to  WNT  pathway  activation  by  WNT3A  stimulation. 
Relative  mRNA  expression  of  Axin2  upon  WNT3A  stimulation  used  as  marker  for  activation  of 
WNT  pathway. 

(C)  DKK1  protein  expression  levels  in  H2087  and  HCC1954-LCC  derivatives  relative  to 
Parental.  Conditioned  media  from  respective  cell  lines  was  assayed  for  DKK1  expression  by 
ELISA.  Data  are  mean  amount  of  secreted  DKK1  protein  ±  S.E.M.  N  =  3  technical  replicates 
per  group. 

(D)  HCC1954-LCC1  cells  (green)  detected  in  the  brain  parenchyma  1  month  post-injection  in 
athymic  nude  mice  express  DKK1  (red). 

(E)  Disseminated  double-positive  DKK1+  (green)  and  human  vimentin+  (red)  cells  were 
detected  in  the  lung  and  kidneys  of  NSG  mice  implanted  with  human  breast  PDX  tumor  chunks 
at  the  mammary  fat  pad.  SOI:  Surgical  Orthotopic  Implantation. 

(F)  Sox2  binding  to  the  DKK1  promoter  assayed  by  Sox2  ChIP  followed  by  qRTPCR  analysis  of 
the  DKK1  promoter  in  H2087  derivatives.  Error  bars  ±  S.D. 

(G)  Dkkl  mRNA  expression  diminished  in  H2087-LCC  cells  depleted  of  Sox2. 

(H)  WNT  pathway  activation  by  WNT3a  stimulation  in  H2087-LCC  cells  depleted  of  DKKl. 
Relative  mRNA  expression  of  Axin2  used  as  marker  for  WNT  pathway  activation. 

(I)  Dye  retention  assay  on  H2087-LCC1  cells  depleted  for  DKK1  assayed  by  flow  cytometry 
after  6  days  of  in  vitro  culture  in  MLM  conditions.  Scale  bars,  50pm  (D  and  E). 


•  Figure  7.  DKK1  enforces  a  quiescent  immune  evasive  state. 

(A)  Growth  of  H2087-LCC  cells  depleted  of  DKK1  in  the  restrictive  environment  of  athymic 
nude  mice.  Data  are  total  number  of  LCC  cells  detected  per  lung  after  3  months  post-tail  vein 
injection  ±  S.E.M.  N  =  5-6  mice  per  group,  scoring  representative  serial  sections  of  the  entire 
lung  for  each  mouse. 

(B)  Relative  mRNA  expression  of  several  NK  cell  activating  ligands  and  death  receptors  in 
H2087-LCC  cells  upon  DKK1  depletion.  Relative  Quantification  is  normalized  to  LCC1  or 
LCC2  Control. 

(C)  Depletion  of  NK  cells  by  anti-asialo-GMl  in  athymic  nude  mice  injected  with  shDKKl 
bearing  H2087-LCC2  cells.  Data  are  whole  body  photon  flux  2  months  post-injection  with 
whisker  plots  representing  minimum  and  maximum  values.  N  =  5-6  mice  per  group.  P  <0.01 
(**),  Mann- Whitney  Test. 

(D)  Growth  of  H2087-LCC  cells  depleted  of  DKK1  in  NSG  mice.  Data  are  whole  body  photon 
flux  2  months  post-injection  with  whisker  plots  representing  minimum  and  maximum  values.  N 
=  5-6  mice  per  group.  P  <0.01  (**),  P  <0.001  (***),  Mann- Whitney  Test. 
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(E)  Representative  immunofluorescence  images  of  metastatic  lesions  in  the  lungs  of  NSG  mice 
in  Figure  7D. 

(F)  Growth  of  H2087-LCC  cells  depleted  of  SOX2  in  NSG  mice.  Data  are  whole  body  photon 
flux  2  months  post-injection  with  whisker  plots  representing  minimum  and  maximum  values.  N 
=  5-6  mice  per  group.  P  <0.05  (*),  P  <0.001  (***),  Mann- Whitney  Test. 

(G)  A  model  summarizing  central  tenets  of  latency  recapitulated  in  our  latency  competent  cancer 
cell  models.  Disseminated  LCC  cells  enriched  for  Sox2/Sox9  stem  cell-like  programs  remain 
associated  with  the  vasculature,  maintain  their  slow-cycling  status  by  inhibition  of  proliferative 
WNT  signaling  through  secretion  of  autocrine  DKK1,  and  avoid  NK  cell  mediated  cell  death. 

Scale  bars,  75pm  (D)  and  50pm  (F). 


•  Figure  8.  Inhibition  of  gap  junction  activity  controls  brain  metastatic  outgrowth. 

(A)  Schematic  summary  of  gap  junction  mediated  anti-dsDNA  response,  production  of  IFNa  and 
TNFa  in  astrocytes,  and  consequent  activation  of  STAT1  and  NF-kB  pathways  in  cancer  cells  to 
support  brain  metastasis. 

(B)  Dye  transfer  from  astrocytes  to  MDA231-BrM2  cells  in  the  presence  of  the  indicated 
concentrations  of  Tonabersat  or  meclofenamate.  n  >  3  independent  experiments. 

(C)  ELISA  of  IFNa  and  TNFa  in  conditioned  media  from  co-cultured  MDA231-BrM2  cell  and 
astrocytes  in  the  presence  of  Tonabersat  (Tona)  or  meclofenamate  (Meclo)  with  indicated 
concentrations.  All  graphs  shown  are  mean  ±  S.E.M.  (n=4  technical  replicates).  n=2  independent 
experiments. 

(D)  Tonabersat  or  meclofenamate  was  administered  daily  starting  one  day  after  cancer  cell 
inoculation  in  mice.  Brain  metastatic  lesions  were  quantified  based  on  BLI.  n=2  independent 
experiments. 

(E)  GFP  staining  of  14-day  brain  metastatic  lesions.  Representative  images  show  large, 
progressive  lesions.  DAPI,  nuclear  staining.  Scale  Bar,  40pm.  n=10  experimental  mice. 

(F)  14  days  after  inoculation  with  MDA231-BrM2  cells  transduced  with  inducible  control,  CX43 
or  PCDH7  shRNAs,  mice  were  treated  with  doxycycline  and  carboplatin,  as  illustrated  in  the 
scheme. 

(G,  H)  Brain  metastatic  lesions  were  quantified  based  on  BLI.  Representative  images  of  matched 
ex  vivo  brain  BLI  and  red  fluorescence  imaging.  n=2  independent  experiments. 

(I,  J)  14  days  after  inoculation  with  MDA231-BrM2  cells,  mice  were  treated  with  Tonabersat, 
meclofenamate,  and  carboplatin.  Following  the  indicated  regimens  (I),  brain  metastatic  lesions 
were  quantified  based  on  BLI  (J).  n=2  independent  experiments. 


•  Publication:  Obenauf,  A.C.,  Zou,  Y.,  Ji,  A.L.,  Vanharanta,  S.,  Shu,  W.,  Shi,  H.,  Kong,  X., 
Bosenberg,  M.C.,  Wiesner,  T.,  Rosen,  N.,  Lo,  R.S.,  and  Massague,  J.  Therapy  induced  tumour 
secretomes  promote  resistance  and  tumour  progression.  Nature  520,  368-372  (2015) 
PMC4507807 


I  Figure:  Figures.pdf 


Figure  1 


Days  Post  Injection 


C 


D  H2087-LCC1 

H2087-LCC1 

H2087-LCC1 

..  HCC1954;LCC1 

CD31  1 

,t  •‘Cpi  y v 

v  c 

ft  N  J 

,  ' 

^  %.  \  7 

O 

Kidney 

Kidney 

Lung 

Brain 

' .  •.*,  .  ^  ** 

■  _ 

%  -•  /  « •  * 

•  •  V  •••'  !  *  '-** 

G 

huVimentin 

-  —  ^  • 

f  •  / 

’%••••’  •  Mt  • _ 

••  *•*’•  *•  *< 

t  T  •  * 1  i  •  1 

4  •  '  *•#  %  - 

•*  •  *  •>  ••  •  _ •  _ 

•••*  ^ 

Kidney 

# 

Kidney  ;.•/  V*  * 

Kidney 

_ _ 

Kidney 

_ 

r.  *•  "*'i 

- .  .  ,■  .  . 

it*  V 

01  r 

V  n  ^  ^ 

U?  V  V* 

*.  -  V  >>  ' 

w  »  '« 

•  *  W  V  * 

«  *  t 

wk  k  w  ‘ 

^ 

?  . 

y 

i  -•.s*.  rc  •  < 

Lung  -V  •  ^ 

Brain  #  *  * 

Lung 

Brain 

_ 

H2087-LCC2 

V.CD31" 


^  >  > 


'Me- 


Kidney 


Figure  2 


A 


B 


C 


EdU  Labeled 
LCC  Cells  l-5*  •*  2  wks 


Quantify 
EdU+ cells  in 
lungs 


Unlabeled  EdU  Labeled 


LCC  Cells 


-|00  |  Parental  |  LCC1 


CO 


D 


E 


F 


H2087-LCC1 

H2087-LCC1 

HCC1954-LCC1 

EdU  •  ,  X 

Ki-67 

•  , 

D14 

* >  ••  _ 

ffh  * 

ft#  ^ 

o 

CD 

O 

$ 

|C 

■Ki67Nes  BKi67pos 

O  75- 
O 
O 
_ i  15-r 

i'l  I  i 

^  o  II  M —  \m 

H2087  HCC1954 


G  H 

—  Parental  — LCC1  — LCC2 


APC  APC 


00 

o 

C\J 

x 


LO 

o> 

O 

O 

X 


Cell  Cycle  States 


Parent 

LCC1 

LCC2 

Parent 


LCC1 


LCC2 


G0/G1  S  G2/M 

1  1  1 

i  r 

i  1 1 

1 1 1 

i  n 

i  rr 

i  i 

i  i 

i  i 

i  i 

i  i 

l l 

- 1 - 1 - 1 - 1 

MRM 

MLM 

MRM 

MLM 

MRM 

MLM 

MRM 

MLM 

MRM 

MLM 

MRM 

MLM 


0  25  50  75  100 


Percentage  of  Cells 


#  of  LCC  Cells/Brain  Relative  mRNA  Levels  PC2: 1 1  %  Variance 


Figure  3 


40- 


20- 


0- 


-20- 


-40- 


Stroma 


CP0 


-60  -40  -20  0  20  40 

PCI :  23%  Variance 
^HumanQMouse  /\  MLM\7mRM 


0  102  1 03  1 04  1 05 


CD44-APC 


PCI :  27%  Variance 

OMouse  /\MLMVMRM 


□  Parental  BLCC1  BLCC2 


H 


=  -r-  (N 

£  O  O 

CO  o  o 

CL  I  I 


Ctrl 


shSox2-1  shSox2-2 


400  i 


00 

o 

C\J 

X 


LO 

ct> 

o 

O 

X 


O) 

C 


300  - 


Sox2 

Sox9  § 

a3 

-r  .  ,  °  200 

Tubulin  q 
O 

_i 

Sox2  o  1 00  - 

Sox9 
Tubulin 


T 


LCC1 


LCC2 


Ctrl 

500  i  shSox9 
400 
300 
200 
100 
0 

HCC1954-LCC1 


Parental 

I 


i  LCC1 

*  t  - 


shSox2-1  shSox2-2 


Parental 


LCC1 


LCC2 


Whole  Body  Photon  Flux 


Figure  4 


A 


B 


C 


H2030-BrM 

H2087-LCC1 

H2087-LCC2 


15  30  45  60 

Days  Post  Injection 


<§  o 

c\j  O 
x  -1 


cvi 

CO  f'v 

°  H 

CM  O 
X  -I 


^1- 

LO 

2o 
o  o 
o  -1 

X 


D17  D30  D60 


D 


E 


F 


G 


|CtrlHPK136|GM1 


■Ctrl  HGM1 


T-nonMet  T-Met 


Death  Receptors  NK  Cell  Activating  Ligands 


Figure  5 


B 


HCC1954-LCC1 
■MRM  MLM 

Relative  mRNA  Levels 
0.0  0.5  1.0  1.5 


CD155 


H2087-LCC1 
■MRM  MLM 

Relative  mRNA  Levels 
0.0  0.5  1.0  1.5 


CO 

~G 


as 

O) 


O) 


> 

o 

< 


CD155 


ULBP1 


ULBP2 


ULBP4 


=5  ULBP5 


05 

i— 

o 

Q_ 


MICA 


ICAM1 


FAS 


£  TRAILR1 


DC 


£  TRAILR2 


C 


D 


_  _  Spleen 

Cancer  cells 

Label  -  e670  Isolate  NK  cells 


) 


MLM 


12  hrs  IL-2 


1:5 

\3  hrs  ^ 


I 

%  Cytolysis 


□  Parental  HLCC1  HLCC2 
60  n 


H2087  HCC1954 


Rel.  DKK1  mRNA  Levels 


Figure  6 


A 


H2087  -  MLM 
Parental  LCC  Cells 
4  2  o  1  2  IBB— 

MYC 
WNT 
NF-kB 
RAS 
TGF-b 
SRC 
ERBB2 
bCAT 
AKT 
E2F1 
PI3K 

MYC 

|  Signature  J  |  mRNA 

z-score  -j  q  -1  z-score 


B 


c 


□  CtrlH  WNT3A 


□  Parental  ■  LCC1  El  LCC2 


1  M  M  0  11  l  "  " 

H2087  HCC1954 


D 


HCC1954-LCC1 

DKK1 

y  ■ 

'U 

‘  Jltf' 

y  i 

Brain 

GFP 

**  6"kki 

F 


SOI:  Breast  PDXs 


D60  -  Harvest  organs 

I 


Section  &  Stain  for 
Disseminated  Tumor 
Cells 


Lung  Kidney 


Vimentin 

DKK1 

r  ^ 

MERGE 

* 

>  \ 

0.30 

0.25 

§0.20 

CD 

£0.15 

_Ctf 

cr  0.10 
0.05 


DKK1  Promoter 
□  Ctrl  |  SOX2 


r 


Pa  LCC1  LCC2 


G 


1.5-. 


shCtrl 

■  shSox2-1  BshSox2-2 


LCC1  LCC2 


□  Ctrl  ■  WNT3A 


shCtrl  shDKK1-2 


Whole  Body  Photon  Flux 


Figure  7 


B 


Ctrl  shDKK1-lHshDKK1-2 


300 

0) 

c 

0 

_l 

1  200 

CD 

O 

o 

O 


100 


o 

=H= 


LCC1 


LCC2 


H2087  ■  LCC1  Ctrl  ■  LCC1  shDKKI 


_c n 
CD 
> 
0 
_l 

< 

z 

dc 


0 

> 

CD 

DC 


LO 

T— 

C\J 

LO 

< 

T— 

Cf) 

T— 

C\J 

LO 

CL 

CL 

CL 

CL 

o 

< 

DC 

DC 

0Q 

CO 

DO 

CQ 

< 

LL 

_l 

_l 

Q 

O 

_l 

3 

_l 

3 

_l 

3 

_l 

3 

o 

< 

DC 

1- 

< 

DC 

1- 

NK  Cell  Activating  Ligands  Death  Receptors 


1 04n 


103- 


102 


Ctrl  |GM1 


x 


shDKKI 
H2087-LCC2 


Ctrl 

shDKKI -1 1  shDKKI -2 


C 

o 

-I — > 

o 

_£Z 

CL 

> 

"0 

O 

DO 


O 


LCC  Cells 

DKK1 

MERGE  b 

.  <57^  a*!.1’*'* 

A]**4 

LCC1 


LCC2 


Ctrl 

shSox2-1 


shSox2-2 


106 


x 

0 


o 

o 


CL  1Q5 

■o 

o 

DO 

0 

O 


104 


0 


4 


LCC1 


LCC2 


*>  Quiescent,  Immune  Evasive  LCC  cell 
*  Proliferating,  Immune  Susceptible  LCC  cell 


%  of  Max 


Figure  8 


B 


Astro -►MDA231-BrM2 


Meclo(|jM):  0  0  0  0  0  0  25  50  100  200  0 

MDA231-BrM2  MDA231-BrM2  +  Astro  Astro 


D 


MDA231 

BrM2 


5  wks  ^ 
Brain  BLI 


1  x106  3 


MDA231 

BrM2 


14  d\ 
Brain  BLI 


F 


I 


H 


P=0.01 


I  Inoculation 

|  |  |  |  |  |  Carboplatin  (i.p.) 

\  Tonabersat  (i.p.) 

\  Tona  +  Carb0 

Meclofenamate  (i.p.) 

I  Mecl°  +  Carb0 

0  14  35  (day) 


Control  Tonabersat  Meclofen. 


LETTER 


doi:10.1038/naturel4336 


Therapy-induced  tumour  secretomes  promote 
resistance  and  tumour  progression 

Anna  C.  Obenauf1,  Yilong  Zou1,2*,  Andrew  L.  Ji1*,  Sakari  Vanharanta1,3,  Weiping  Shu1,  Hubing  Shi4,  Xiangju  Kong4, 
Marcus  C.  Bosenberg5,6,  Thomas  Wiesner7,  Neal  Rosen8,  Roger  S.  Lo4  &  Joan  Massague1 


Drug  resistance  invariably  limits  the  clinical  efficacy  of  targeted 
therapy  with  kinase  inhibitors  against  cancer1,2.  Here  we  show  that 
targeted  therapy  with  BRAF,  ALK  or  EGFR  kinase  inhibitors  induces 
a  complex  network  of  secreted  signals  in  drug-stressed  human  and 
mouse  melanoma  and  human  lung  adenocarcinoma  cells.  This 
therapy-induced  secretome  stimulates  the  outgrowth,  dissemination 
and  metastasis  of  drug-resistant  cancer  cell  clones  and  supports  the 
survival  of  drug-sensitive  cancer  cells,  contributing  to  incomplete 
tumour  regression.  The  tumour-promoting  secretome  of  melanoma 
cells  treated  with  the  kinase  inhibitor  vemurafenib  is  driven  by  down- 
regulation  of  the  transcription  factor  FRA1.  In  situ  transcriptome 
analysis  of  drug-resistant  melanoma  cells  responding  to  the  regres¬ 
sing  tumour  microenvironment  revealed  hyperactivation  of  several 
signalling  pathways,  most  prominently  the  AKT  pathway.  Dual  inhi¬ 
bition  of  RAF  and  the  PI(3)K/AKT/mTOR  intracellular  signalling 
pathways  blunted  the  outgrowth  of  the  drug-resistant  cell  popula¬ 
tion  in  BRAF  mutant  human  melanoma,  suggesting  this  combina¬ 
tion  therapy  as  a  strategy  against  tumour  relapse.  Thus,  therapeutic 
inhibition  of  oncogenic  drivers  induces  vast  secretome  changes  in 
drug-sensitive  cancer  cells,  paradoxically  establishing  a  tumour  micro¬ 
environment  that  supports  the  expansion  of  drug-resistant  clones, 
but  is  susceptible  to  combination  therapy. 

Kinase  inhibitors  such  as  vemurafenib,  erlotinib  or  crizotinib  have 
shown  clinical  efficacy  in  melanoma  with  BRAF  mutations,  or  in  lung 
adenocarcinoma  with  EGFR  mutations  or  ALK  translocations,  respec¬ 
tively3'6.  Although  complete  responses  are  rare,  the  vast  majority  of 
patients  show  partial  tumour  regression  or  disease  stabilization.  How¬ 
ever,  drug  resistance  invariably  develops,  and  most  patients  progress 
within  6-12  months3'16,  representing  a  common  complication  of  tar¬ 
geted  therapies  that  hampers  long-term  treatment  success.  The  rapid 
emergence  of  clinical  drug  resistance  maybe  facilitated  by  a  small  num¬ 
ber  of  pre-existing  cancer  cells  that  are  intrinsically  resistant  or  poised 
to  adapt  to  drug  treatment  quickly17'19.  How  these  minority  clones  of 
drug- resistant  cells  react  to  the  marked  changes  in  the  microenviron¬ 
ment  during  tumour  regression  is  not  known.  A  better  understanding 
of  this  process  could  lead  to  treatments  that  improve  the  efficacy  of  cur¬ 
rent  targeted  anti-cancer  drugs. 

To  model  therapeutic  targeting  of  heterogeneous  tumour  cell  popu¬ 
lations  in  vivo ,  we  mixed  a  small  percentage  of  vemurafenib -resistant 
A3 75  human  melanoma  cells  (A375R),  labelled  with  a  TK-GFP-luciferase 
(TGL)  vector,  together  with  mostly  non-labelled,  vemurafenib -sensitive 
A3 75  cells,  and  injected  the  admixture  (A375/A375R,  99.95/0.05%)  sub¬ 
cutaneously  in  mice  (Extended  Data  Fig.  la).  After  the  tumours  were 
established,  we  treated  the  mice  with  vemurafenib  or  vehicle,  and  mon¬ 
itored  the  growth  of  resistant  cells  by  bioluminescent  imaging  (BLI) 
in  vivo  (Fig.  la).  Although  vemurafenib  treatment  decreased  the  vol¬ 
ume  of  sensitive  tumours  (A375  alone)  (Extended  Data  Fig.  lb),  the 


number  of  admixed  resistant  cells  in  regressing  tumours  (A375/A375R) 
significantly  increased  compared  to  vehicle-treated  controls  (Fig.  lb). 
Green  fluorescent  protein  (GFP)  staining  confirmed  increased  num¬ 
bers  of  resistant  cells  in  regressing  tumours,  and  EdU  or  BrdU  staining 
confirmed  their  increased  proliferation  rate  compared  to  the  vehicle- 
treated  controls  (Fig.  lc  and  Extended  Data  Fig.  lc,  d).  Tumours  com¬ 
prising  only  resistant  cells  showed  no  growth  difference  when  treated 
with  vehicle  or  vemurafenib  (Fig.  Id),  indicating  that  the  growth  advan¬ 
tage  of  resistant  cells  in  regressing  tumours  was  not  caused  by  direct 
effects  of  vemurafenib  on  cancer  or  stromal  cells. 

Treatment  of  mixed  A3 75  and  A375R  tumours  with  dabrafenib,  an¬ 
other  BRAF  inhibitor  (RAFi),  or  doxycycline-induced  knockdown  of 
BRAF  had  similar  effects  (Extended  Data  Fig.  le-g).  In  line  with  these 
findings,  A375R  cells  co-implanted  with  other  vemurafenib -sensitive 
melanoma  cell  lines  (C0I08OO,  LOX  and  UACC62)  also  showed  an  up 
to  eightfold  growth  increase  compared  to  vehicle- treated  control  groups 
(Fig.  le).  Growth  acceleration  of  the  resistant  population  in  a  regressing 
tumour  was  also  observed  in  the  patient- derived8  melanoma  cell  line 
M249  and  its  vemurafenib-resistant  derivative  M249R4,  driven  by  an 
NRAS  mutation,  a  clinically  relevant  resistance  mechanism  (Fig.  le  and 
Extended  Data  Fig.  lh).  In  immunocompetent  mice,  vemurafenib  treat¬ 
ment  of  tumours  formed  by  melanoma  cell  lines  derived  from 
Braf600E  Cdkn2a~/~  Pten mice  (YUMM1.1,  YUMM1.7)  also  pro¬ 
moted  growth  of  the  admixed  vemurafenib-resistant  cells  (YUMM1 .7R, 
B16)  (Extended  Data  Fig.  li,  j). 

Crizotinib-  or  erlotinib -treated  mice  containing  tumours  formed  by 
AEK-driven  (H3122)  or  EGER-driven  (HCC827)  human  lung  adeno¬ 
carcinoma  cells,  respectively,  admixed  with  minority  clones  of  intrins¬ 
ically  resistant  cells  from  the  same  cell  lineage  (lung  adenocarcinoma 
cells  H2030  and  PC9)  or  melanoma  cells  (A375R)  also  led  to  increased 
outgrowth  of  the  resistant  cells  (Fig.  le  and  Extended  Data  Fig.  lk-m). 
Local  growth  acceleration  of  resistant  cells  in  the  regressing  subcutan¬ 
eous  tumours  resulted  in  higher  lung  metastatic  burden  (Fig.  If).  Thus, 
drug- resistant  cancer  cells  benefit  from  therapeutic  targeting  of  sur¬ 
rounding  drug- sensitive  cells. 

Circulating  tumour  cells  can  infiltrate  and  colonize  tumours.  This 
phenomenon,  termed  self- seeding20,  may  contribute  to  the  distribution 
of  resistant  clones  to  several  metastatic  sites.  Mice  implanted  with  sen¬ 
sitive  A375  tumours  were  treated  with  vehicle  or  vemurafenib,  and  in- 
tracardially  injected  with  TGL-labelled  A375R  cells  (Fig.  lg).  A375R  cells 
were  more  efficiently  attracted  to  vemurafenib-treated  regressing  tumours 
compared  to  vehicle- treated  controls,  with  95%  (2 1  out  of  22)  and  12.5% 
(2  out  of  16)  efficiency,  respectively,  exhibiting  substantial  accumulation 
of  resistant  cells  in  regressing  tumours  by  day  5  (Fig.  lg  and  Extended 
Data  Fig.  In).  To  evaluate  the  contribution  of  seeding  by  resistant  cir¬ 
culating  tumour  cells  to  disease  relapse,  we  intracar dially  injected  re¬ 
sistant  A375r  cells  or  vehicle  into  tumour-bearing  mice  and  compared 
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Figure  1  |  The  regressing  tumour 
microenvironment  stimulates  the  outgrowth, 
infiltration  and  metastasis  of  drug-resistant 
clones,  a,  Schematic  of  the  experimental  set-up. 
b,  Bioluminescent  signal  of  drug- resistant  A375R- 
TGL  cells  in  vemurafenib-sensitive,  A375  tumours, 
treated  with  vehicle  or  vemurafenib  for  5  days 
(vehicle,  n  —  36;  vemurafenib,  n  =  15  tumours). 

D,  day.  c,  EdU  incorporation  in  A375R-TGL  cells  in 
A375/A375r-TGL  tumours  treated  with  vehicle  or 
vemurafenib  for  4  days,  as  determined  by  FACS 
(vehicle,  n  —  8;  vemurafenib,  n  —  6  tumours), 
d,  Bioluminescent  signal  of  A375R-TGL  tumours 
alone,  treated  with  vehicle  or  vemurafenib  for 
5  days  (vehicle,  n  —  38;  vemurafenib,  n—  15 
tumours),  e,  Bioluminescent  signal  of  TGL- 
expressing  drug-resistant  cancer  cells  (A375R, 
M249R4,  PC9  and  H2030)  in  drug- sensitive 
tumours  (C0I08OO,  LOX,  UACC62,  M249,  H3122 
and  HCC827)  treated  with  vehicle  or  drugs 
(vemurafenib,  crizotinib  and  erlotinib)  for  5  days 
( n  (from  left  to  right  on  the  graph)  =  6, 7, 12, 12, 9, 
9,  25,  26,  9,  12,  12,  12,  16  and  11  tumours), 
f,  Spontaneous  lung  metastasis  by  A375R  cells  in 
mice  bearing  A375/A375R-TGL  tumours  treated 
with  vehicle  or  vemurafenib  (10  days),  visualized 
by  BLI  ( n  —  4).  g,  Seeding  of  A375R-TGL  cells  from 
the  circulation  to  unlabelled,  subcutaneous  (SC) 
A3 75  tumours  of  mice  treated  with  vehicle  or 
vemurafenib.  Signal  in  the  tumour  was  quantified 
by  BLI  (vehicle,  n  —  30;  vemurafenib,  n  =  34 
tumours;  three  independent  experiments 
combined).  IC,  intracardiac,  h,  Treatment 
response,  determined  by  tumour  size,  of 
subcutaneous  A375  tumours  allowed  to  be  seeded 
by  A375r-TGL  cells  from  the  circulation  or  mock 
injected  (vehicle,  n  —  16;  vemurafenib,  n  —  8 
tumours).  Data  in  b-e,  g,  h  are  mean  and  s.e.m; 
in  f  the  centre  line  is  median,  whiskers  are 
minimum  and  maximum  values.  *P  <  0.05, 

**P  <  0.01,  ***P  <  0.001,  two-tailed  Mann- 
Whitney  U  test.  NS,  not  significant. 


the  tumour  volume  during  vemurafenib  treatment  (Fig.  lh).  Whereas 
the  unseeded  tumours  in  the  control  group  showed  extensive  tumour 
regression,  seeding  by  A3 75R  cells  led  to  rapid  tumour  relapse  (Fig.  lh). 
These  results  suggest  that  tumours  regressing  on  targeted  therapy  are 
potent  attractors  of  resistant  circulating  tumour  cells  that  may  contri¬ 
bute  to  rapid  tumour  progression. 

Tumours  consist  of  a  complex  microenvironment  composed  of  im¬ 
mune,  stromal  and  cancer  cells21.  Soluble  mediators  from  this  micro¬ 
environment  can  foster  cancer  growth  and  therapy  resistance13,14,22'24. 
Considering  that  drug- sensitive  cancer  cells  are  the  main  population 
affected  by  targeted  therapy,  we  proposed  that  signals  derived  from  sen¬ 
sitive  cancer  cells  in  response  to  kinase  inhibitors  drive  the  outgrowth 
of  drug-resistant  cells.  To  test  this  hypothesis,  we  established  an  in  vitro 
co-culture  system  and  monitored  the  growth  of  TGL- expressing  res¬ 
istant  cells  (A375r,  H2030)  in  the  absence  or  presence  of  sensitive  cells 
treated  with  kinase  inhibitors  or  vehicle  (Fig.  2a).  Mimicking  our 
in  vivo  findings,  co-culture  with  vemurafenib-,  crizotinib-  or  erlotinib- 
treated  sensitive  cells  significantly  enhanced  the  growth  of  resistant  can¬ 
cer  cells  (Fig.  2a  and  Extended  Data  Fig.  2a-c). 

We  derived  conditioned  media  (CM)  from  vemurafenib-sensitive 
melanoma  cells  cultured  in  the  absence  (CM-vehicle)  or  presence  of 
vemurafenib  (CM -vemurafenib).  CM -vemurafenib  accelerated  the  prolif¬ 
eration  of  drug-resistant  cells,  with  different  clinically  relevant  resistance 
mechanisms,  as  determined  by  cell  viability  assays  and  Ki67  staining 
(Fig.  2b  and  Extended  Data  Fig.  2d-f).  Similarly,  conditioned  media 
from  crizotinib-  or  erlotinib -treated  sensitive  lung  adenocarcinoma  cells 
stimulated  proliferation  of  lung  adenocarcinoma  cells  with  intrinsic  or 
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acquired  resistance  (Fig.  2c)  and  across  different  cell  lineages  (Extended 
Data  Fig.  2g).  In  addition,  CM -vemurafenib  elicited  increased  cell  mi¬ 
gration  in  transwell  migration  and  monolayer  gap-closing  assays  (Fig.  2d 
and  Extended  Data  Fig.  2h-k).  CM- vemurafenib  was  also  active  on 
vemurafenib-sensitive  cancer  cells,  increasing  survival  and  suppressing 
the  apoptotic  caspase  activity  up  to  100-fold  in  these  cells  when  treated 
with  vemurafenib  in  vitro  (Fig.  2e,  f).  Because  all  biologically  active  con¬ 
ditioned  media  was  collected  before  cell  death  or  senescence,  it  is  likely 
that  the  secretome  is  actively  produced  as  a  result  of  oncogene  inhibi¬ 
tion  (Extended  Data  Fig.  21,  m).  These  results  demonstrate  that  BRAF , 
ALK  and  EGFR  mutant  cells  respond  to  therapeutic  stress  under  tar¬ 
geted  therapy  by  secreting  factors  that  support  the  survival  of  drug- 
sensitive  cells  and  accelerate  the  growth  of  drug- resistant  minority  clones. 
The  effects  of  this  reactive  secretome  may  augment  previously  reported 
resistance  mechanisms  including  relief  of  feedback  inhibition  of  intra¬ 
cellular  signalling11,25,  upregulation  of  receptor  tyrosine  kinases26,  or 
the  supply  of  stromal  cytokines14  that  protect  the  drug- sensitive  cells. 

To  identify  relevant  components  and  regulators  of  the  reactive  se¬ 
cretome,  we  analysed  gene  expression  changes  in  sensitive  A3 75  mel¬ 
anoma  cells  at  different  time  points  after  vemurafenib  exposure  in  vitro. 
After  6h  on  vemurafenib,  473  genes  showed  altered  expression,  and 
pathway  analysis  revealed  that  these  genes  were  enriched  for  transcrip¬ 
tional  regulators  (Fig.  3a,  b,  Extended  Data  Fig.  3a,  b  and  Supplemen¬ 
tary  Table  1).  After  48  h,  more  than  one- third  of  the  transcrip  tome  was 
differentially  expressed  (>5,000  genes;  405  genes  encoding  for  proteins 
in  the  extracellular  region,  Gene  Ontology  (GO)  accession  0005576), 
significantly  overlapping  with  the  gene  expression  changes  of  A3 75 
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Figure  2  |  The  secretome  of  RAF  and  ALK  inhibitor-treated  tumour  cells 
increases  proliferation  and  migration  of  drug-resistant  cells  and  supports 
the  survival  of  drug-sensitive  cells,  a,  Schematic  (left)  and  representative 
BLI  images  (right)  after  7  days  of  co-culture.  Average  fold  change  (FC)  of  BLI 
signal  from  A375R-TGL  cells  in  vemurafenib -treated  wells  relative  to  vehicle- 
treated  control  wells  is  depicted  on  the  right  {n  =  4  biological  replicates), 
b,  c,  Conditioned  media  (CM)  was  derived  from  drug-sensitive  cells,  treated 
with  vehicle,  vemurafenib  or  crizotinib.  Drug-resistant  cells  were  grown  in 
this  conditioned  media  and  the  cell  number  was  determined  on  day  3. 

Drug- sensitive  and  drug-resistant  cell  lines  and  drugs  used  to  generate 


conditioned  media  as  indicated,  n  =  3  (b)  and  6  (c)  biological  replicates, 
d,  Schematic  diagram  of  the  migration  assay  (top)  and  relative  migration  of 
A375r  cells  towards  conditioned  media  from  different  sources  as  indicated 
(bottom,  n  =  10  fields  of  vision  (FOV)).  ****p  <  0.0001,  two-tailed  Mann- 
Whitney  U  test,  e,  Survival  assay  of  drug-sensitive  A375  cells  cultured  in 
conditioned  media  and  treated  with  vemurafenib,  assessed  on  day  3  (n  =  3 
biological  replicates),  f,  Apoptosis  rate  of  A375  cells  cultured  in  conditioned 
media  and  treated  with  vemurafenib  (3  pM)  {n  —  3  biological  replicates). 
Data  are  mean  and  s.e.m. 
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Figure  3  |  FRA1  downregulation  during  RAFi  treatment  drives  the  reactive 
secretome.  a,  Principal  component  (PC)  analysis  of  drug- sensitive  A375 
cells  treated  in  vitro  with  vehicle  or  vemurafenib  for  6  or  48  h.  b,  Volcano 
plots  show  genes  significantly  deregulated  by  vemurafenib  treatment  after  6  h 
(left)  or  48  h  (right).  Transcription  factors  (TF)  and  gene  products  in  the 
extracellular  region  are  depicted  in  green  (downregulated)  and  red 
(upregulated)  {n  —  3  tumours).  ^adj>  adjusted  P  value,  c,  Relative  mRNA 
levels  of  FRA1  during  vemurafenib  exposure  [0.1-1  pM].  d,  Representative 
immunofluorescence  staining  of  A375/A375R  tumours  for  GFP  (A375R,  green) 
and  FRA1  (red)  after  vehicle  or  vemurafenib  treatment  (5  days).  DAPI, 
4\6-diamidino-2-phenylindole.  Scale  bars,  50  pm.  e,  Top,  representative 


immunofluorescence  staining  for  FRA1  (red)  of  melanoma  biopsy  sections  of 
patient  1.  Original  magnification,  X20.  Bottom,  nuclear  FRA1  staining  was 
quantified  in  three  melanoma  patients  before  (B)  and  early-on  therapy.  RAFi 
and  MEKi  denote  RAF  and  MEK  inhibitors,  respectively,  f,  Bioluminescent 
signal  of  A375R-TGL  cells  6  days  after  subcutaneous  co-implantation  with 
A375  cells  expressing  control  (shCtrl)  or  two  independent  FRA1  shRNAs 
(shFRAl-1  andshFRAl-2)  ( n  —  16  tumours),  g,  Seeding  of  A375R-TGL  cells  to 
unlabelled  tumours  expressing  control  or  two  independent  shRNAs  for  FRA1 , 
determined  by  BLI  (vehicle,  n  =  10;  shFRAl-1,  n  =  10;  shFRAl-2,  n  —  8 
tumours).  Data  are  mean  and  s.e.m.  *P  <  0.05,  **P  <  0.01,  ****p  <  0.0001, 
Student’s  t-test. 
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tumours  in  vivo  after  5  days  of  vemurafenib  treatment  (Fig.  3a,  b  and 
Extended  Data  Fig.  3c).  Similar  extensive  gene  expression  changes 
were  observed  in  C0I08OO  and  UACC62  melanoma  cells  treated  with 
vemurafenib  and  H3122  lung  adenocarcinoma  cells  treated  with  cri- 
zotinib  (Extended  Data  Fig.  3d).  Despite  different  cell  lineages,  differ¬ 
ent  oncogenic  drivers,  and  different  targeted  therapies  we  observed  a 
significant  overlap  between  the  secretome  of  melanoma  and  lung  ade¬ 
nocarcinoma  cells  (P<9.11  X  10-5)  (Extended  Data  Fig.  3e-h  and 
Supplementary  Table  1).  Furthermore,  changes  in  the  secretome  of 
vemurafenib -sensitive  melanoma  cells  coincided  with  changes  in  the 
immune  cell  composition  (Extended  Data  Fig.  4a,  b),  and  with  changes 
of  soluble  mediators  derived  from  murine  stromal  cells  such  as  IGF1 
and  HGF  (Extended  Data  Fig.  4c,  d).  These  data  indicate  a  therapy- 
induced  secretome  (TIS),  a  response  that  consists  of  many  up-  and 
downregulated  secreted  factors,  permeates  the  regressing  tumour  mi¬ 
croenvironment  and  stimulates  cancer  cells,  probably  also  stromal  cells. 

To  identify  molecular  drivers  of  the  A375-TIS  in  response  to  vemur¬ 
afenib,  we  integrated  the  data  of  differentially  expressed  transcription 
factors  after  6  h  of  vemurafenib  treatment  with  the  transcription  factor 
binding  motifs  that  were  enriched  at  the  promoters  of  differentially 
expressed  genes  in  the  secretome  after  48  h  (Fig.  3a,  b).  This  analysis 
highlighted  FRA1  (also  known  FOSL1),  a  member  of  the  API  tran¬ 
scription  factor  complex  and  effector  of  the  ERK  pathway27,  as  one  of 
the  putative  upstream  regulators  of  the  TIS  (Extended  Data  Fig.  5a). 
FRA1  was  downregulated  in  all  drug- sensitive  cells,  but  not  in  resistant 
cells,  treated  with  vemurafenib,  crizotinib  and  erlotinib  (Fig.  3c,  d  and 
Extended  Data  Fig.  5b-d).  Biopsies  from  melanoma  patients  early  during 
RAFi  treatment  confirmed  RAFi-induced  FRA1  downregulation  in  clin¬ 
ical  samples  (Fig.  3e,  Extended  Data  Fig.  5e  and  Extended  Data  Table  1). 

To  test  the  functional  role  of  FRA  1  in  modulating  the  TIS,  we  used 
RNA  interference  (RNAi)  to  inhibit  FRA1  expression.  Co-culture  and 
conditioned  media  assays  using  A3 75  cells  expressing  short  hairpin 


RNAs  targeting  FRA1  (shFRAl)  showed  similar  growth-accelerating  and 
chemotactic  activity  on  A375R  cells  as  vemurafenib  treatment  (Extended 
Data  Fig.  6a-d).  In  line  with  these  results,  FRA1  knockdown  in  A375 
cells  induced  transcriptional  changes  similar  to  those  induced  by  vemur¬ 
afenib  (Extended  Data  Fig.  6e).  A375R  cells  co -implanted  with  A3 75  or 
UACC62  cells  expressing  shFRAl  also  demonstrated  increased  growth 
in  vivo  (Fig.  3f  and  Extended  Data  Fig.  6f).  A375-shFRAl  tumours 
attracted  significantly  more  resistant  cells  from  the  circulation  than 
tumours  expressing  the  control  vector  (Fig.  3g).  Thus,  FRA1  down- 
regulation  drives  the  induction  of  the  tumour-promoting  secretome  of 
vemurafenib -treated  cancer  cells. 

To  determine  the  effect  of  the  reactive  secretome  on  the  drug- 
resistant  tumour  subpopulation  in  a  regressing  tumour,  we  expressed 
the  ribosomal  protein  LlOa  (RPLIOa)  fused  to  enhanced  green  fluor¬ 
escent  protein  (eGFP-RPLIOa)  in  A375R  cells,  allowing  the  specific 
retrieval  of  transcripts  from  A375R  cells  by  polysome  immunopreci- 
pitation  for  subsequent  RNA- sequencing  (RNA-seq)  analysis28  (Fig.  4a). 
In  line  with  the  in  vivo  phenotype  of  accelerated  growth,  the  gene  ex¬ 
pression  pattern  of  resistant  cells  in  the  regressing  microenvironment 
was  enriched  for  biological  processes  involved  in  cell  viability,  prolif¬ 
eration  and  cell  movement  (Extended  Data  Fig.  7a).  Pathway  analysis 
of  the  expression  data  suggested  activation  of  several  pathways  includ¬ 
ing  PI(3)K/AKT,  BMP-SMAD  and  NFkB  (Fig.  4b).  The  hyperactivity 
of  the  PI(3)K/AKT  pathway  in  this  context  also  suggested  a  potential 
vulnerability  of  the  cells  to  PI(3)K/mTOR  inhibitors  (Extended  Data 
Fig.  7b).  The  pathway- analysis-based  prediction  ofPI(3)K/AKT  activa¬ 
tion  was  also  reflected  at  the  protein  level  in  both  resistant  and  sensitive 
cells  in  the  presence  of  CM -vemurafenib  in  vitro  and  under  vemur¬ 
afenib  treatment  in  vivo  (Fig.  4c  and  Extended  Data  Fig.  7c,  d).  More¬ 
over,  PI(3)K/AKT  emerged  as  the  dominant  TIS  responsive  pathway 
in  a  targeted  immunoblot  analysis  of  survival  pathways  in  vitro  (Ex¬ 
tended  Data  Fig.  7e). 
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Figure  4  |  The  therapy-induced  secretome  in  melanoma  promotes  relapse 
by  activating  the  AKT  pathway  in  resistant  cells,  a,  Schematic  diagram 
showing  the  isolation  of  polysome- associated  transcripts  from  resistant  cells  by 
translating  ribosome  affinity  profiling  (TRAP)  from  tumours  during  treatment. 
IP,  immunoprecipitation.  b,  Ingenuity  upstream  regulator  analysis  of  gene 
expression  profiles  from  A375R  cells  responding  to  a  regressing  tumour 
microenvironment  (5  days  of  treatment;  n  —  3  tumours),  c,  Phosphorylation 
status  of  AKTS473  (pAKT)  in  A375R  cells,  stimulated  for  15  min  with  various 
conditioned  media,  as  indicated  by  immunoblotting.  tAKT,  total  AKT. 
d,  Phosphorylation  status  of  AKTS473  in  A375R  cells  after  stimulation  with 


positive  regulators  of  the  AKT  pathway,  upregulated  in  the  melanoma  TIS; 
ANGPTL7  (5  jig  ml-1,  30  min;  upregulated  in  A375,  C0I08OO,  UACC62), 
PDGFD  (lOngml-1,  10  min;  upregulated  in  C0I08OO),  EGF  (lOngml-1, 

10 min;  upregulated  in  A375)  and  IGF1  (lOngml-1,  10 min;  upregulated  in 
UACC62).  e,  Mice  bearing  A375/A375R-TGL  tumours  were  treated  with 
drugs,  and  growth  of  A375R  cells  was  followed  by  BLI  (vehicle,  n—  14; 
vemurafenib,  n  —  16;  vemurafenib  and  BEZ235,  n  —  16;  vemurafenib  and 
MK2206,  n  =  8  tumours),  f,  Graphical  summary  of  the  findings.  Data  are  mean 
and  s.e.m.  P  values  calculated  using  a  two-tailed  Mann-Whitney  U  test. 
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The  TIS  contained  many  mediators  directly  or  indirectly  activating 
the  AKT  pathway.  Positive  mediators  that  were  upregulated  during 
therapy  included  IGF1,  EGF,  ANGPTL7  and  PDGFD,  each  of  which 
activated  the  AKT  pathway  in  vitro  (Fig.  4d).  IGF1,  one  of  the  most 
potent  activators  of  the  AKT  pathway,  is  also  abundantly  expressed  in 
the  tumour  stroma  and  is  further  upregulated  during  targeted  therapy 
(Extended  Data  Figs  4c  and  7 f).  In  addition,  levels  of  IGFBP3,  a  nega¬ 
tive  regulator  of  IGF  1,  were  markedly  reduced  in  the  TIS  of  all  inves¬ 
tigated  cell  lines,  favouring  increased  AKT  pathway  activation  in  the 
presence  of  IGF1  and  stimulation  of  proliferation  of  resistant  cells 
in  vivo  (Extended  Data  Fig.  7f-k). 

To  test  the  role  of  AKT  activation  as  a  mediator  of  TIS-induced  tu¬ 
mour  proliferation,  we  combined  vemurafenib  with  AKT/PI(3)K/mTOR 
inhibitors.  In  co- culture  and  proliferation  experiments  using  condi¬ 
tioned  media,  dual  inhibition  of  the  MAPK  and  AKT  pathway  dimin¬ 
ished  the  growth  benefit  of  the  TIS  (Extended  Data  Fig.  8a,  b).  We  then 
treated  mice  with  A375/A375R  or  A375R  tumours  with  vemurafenib 
and  AKT  (MK2206)  or  PI(3)K/mTOR  inhibitors  (BEZ235).  The  com¬ 
bined  inhibition  of  MAPK  and  PI(3)K/AKT/mTOR  pathways  signifi¬ 
cantly  blunted  the  outgrowth  of  vemurafenib -resistant  cells  in  the  A3 75/ 
A375r  tumours  (Fig.  4e).  The  growth  inhibition  was  specific  for  the 
amplified  proliferation  in  the  regressing  tumour  microenvironment  and 
had  no  effects  on  the  growth  of  resistant  cells  alone  (Extended  Data 
Fig.  8c).  Furthermore,  the  outgrowth  of  resistant  A375R  cells  in  tumour 
seeding  assays  was  significantly  reduced  when  regressing  tumours 
were  co-treated  with  BEZ235  (Extended  Data  Fig.  8d).  Thus,  the 
TIS-induced  proliferation  is  susceptible  to  therapeutic  targeting. 

The  limited  effectiveness  of  targeted  therapies  has  been  attributed  to 
intracellular  feedback  loops  and  specific  cytokines  that  support  the 
survival  of  drug- sensitive  cells.  From  these  residual  tumours,  clones 
emerge  that  are  intrinsically  resistant  to  targeted  therapy  and  are  ulti¬ 
mately  responsible  for  clinical  relapse.  Our  work  demonstrates  that 
targeted  inhibition  of  a  cancer  driver  pathway  can  paradoxically  pro¬ 
mote  these  two  aspects  of  drug  resistance  via  induction  of  a  complex, 
reactive  secretome.  This  TIS  not  only  enhances  the  survival  of  drug- 
sensitive  cells,  but  also  acutely  accelerates  the  expansion  and  dissemi¬ 
nation  of  drug-resistant  clones.  Rather  than  a  cell  death  by-product29,30, 
the  TIS  is  a  live-cell  response  to  inhibition  of  an  oncogenic  driver  path¬ 
way,  mediated  by  a  concrete  transcriptional  program,  and  defined  by 
specific  alterations  of  intracellular  signalling  networks  (Fig.  4f). 

Our  identification  of  AKT  signalling  as  a  mediator  of  TIS-induced 
tumour  progression  in  BRAF- driven  melanoma  is  in  line  with  AKT  acti¬ 
vation  in  tumours  observed  in  the  clinic  during  vemurafenib  treatment16. 
Patients  treated  with  BRAF  inhibitor  rarely  show  full  tumour  regres¬ 
sion3,4,  and  the  remaining  drug-responsive  tumour  cells  may  remain  a 
source  of  TIS  for  the  duration  of  the  treatment.  Our  results  provide  a 
rationale  for  combining  PI(3)K/AKT/mTOR  pathway  inhibitors  with 
inhibitors  of  the  MAPK  pathway  in  the  treatment  of  these  tumours. 
However,  the  breadth  of  the  TIS  and  the  generality  of  our  findings 
across  different  cell  lineages,  drugs  (vemurafenib,  crizotinib  and  erlo- 
tinib),  and  resistance  mechanisms  suggest  that  durable  responses  may 
require  the  combination  of  this  type  of  agents  with  a  radically  different 
therapeutic  modality. 

Online  Content  Methods,  along  with  any  additional  Extended  Data  display  items 
and  Source  Data,  are  available  in  the  online  version  of  the  paper;  references  unique 
to  these  sections  appear  only  in  the  online  paper. 
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METHODS 

Cell  culture.  A375,  M249  (ref.  8)  and  B16  cells  were  cultured  in  DMEM  media; 
C0I08OO,  UACC62,  SKMEL239-clone3,  LOX,  PC9,  H2030,  H3122  and  HCC827 
cells  were  cultured  in  RPMI  media.  YUMM1.1  and  YUMM1.7  were  cultured  in 
DMEM/F12  media.  GPG29  and  293T  cells  were  used  for  retrovirus  and  lentivirus 
production,  respectively.  Both  were  maintained  in  DMEM  media.  All  media  contained 
10%  FBS,  2  mM  L-glutamine,  100 IU  ml-1  penicillin/streptomycin  and  1  pgml-1 
amphotericin  B,  the  media  for  GPG29  contained  in  addition  0.3  mg  ml-1  G418, 
20  ng  ml  - 1  doxycycline  and  2  pg  ml  - 1  puromycin.  All  cells  were  grown  in  a  humid¬ 
ified  incubator  at  37  °C  with  5%  C02  and  were  tested  regularly  for  mycoplasma 
contamination.  All  cell  lines  used  were  negative  for  mycoplasma. 

To  generate  vemurafenib-resistant  melanoma  cell  lines,  vemurafenib- sensitive 
cell  lines  were  seeded  at  low  density  and  exposed  to  1  -3  pM  vemurafenib  (LC-Labs). 
After  approximately  8  weeks  of  continuous  vemurafenib  exposure,  we  derived  re¬ 
sistant  cell  clones  that  were  maintained  on  vemurafenib  ( 1  pM  vemurafenib  for 
M249R4,  Co1o800r,  LOXr,  UACC62r;  2  pM  vemurafenib  for  A375R,  YUMM1.7R). 
The  same  protocol  was  performed  to  generate  a  crizotinib-resistant  cell  line  from 
H3 122  lung  adenocarcinoma  cells,  which  were  selected  and  maintained  with  300  nM 
crizotinib.  Drug-sensitive  and  resistant  melanoma  cell  lines  from  A375,  C0I08OO, 
UACC62  and  YUMM1.7  and  the  drug  sensitive  lung  adenocarcionoma  cell  lines 
H3122  and  HCC827  were  exposed  to  increasing  doses  of  vemurafenib  and  the 
number  of  cells  was  determined  after  3  days  and  pERK  levels  after  1  h  of  vemur¬ 
afenib,  crizotinib  or  erlotinib  exposure  (Extended  Data  Fig.  9a-j).  Receptor  status 
was  determined  by  western  blot  and  showed  an  increase  in  EGFR  expression  levels 
in  all  resistant  lines  examined  as  well  as  an  increase  in  MET  receptor  expression  in 
A375r  and  UACC62R  cells  compared  to  their  parental,  drug- sensitive  cells  (Ex¬ 
tended  Data  Fig.  9k). 

For  co-culture  assays  sensitive  cells  were  plated  in  12-well  or  24-well  plates  and 
allowed  to  adhere  overnight  in  regular  growth  media.  Media  was  then  replaced  with 
low  serum  (2%  FBS)  media  containing  vehicle,  0.1  pM  vemurafenib,  0.3  pM  cri¬ 
zotinib  or  0.01  pM  erlotinib.  For  control  wells  media  containing  vehicle  or  0.1  pM 
vemurafenib,  0.3  pM  crizotinib,  or  0.01  pM  erlotinib  was  plated  at  the  same  time. 
After  48  h,  TGL-expressing,  resistant  cells  were  plated  on  top  of  the  vehicle/ drug 
treated  cells  or  in  media-only  control  wells.  Media  containing  vehicle/drug  was  re¬ 
plenished  every  48 h.  After  7 days,  luciferin  [150  pgml-1]  was  added  to  the  wells 
and  luciferase- signal  of  resistant  cells  was  determined  by  BLI  using  a  Xenogen  Spec¬ 
trum  imaging  machine  (Perkin  Elmer) .  Co-culture  experiments  were  independently 
performed  at  least  twice  and  a  representative  experiment  is  shown. 

To  generate  conditioned  media,  2.3  X  106  and  6.4  X  106  drug-sensitive  cells  were 
plated  on  15-cm  dishes  in  regular  growth  media  and  allowed  to  adhere  overnight. 
The  media  was  then  replaced  by  low  serum  media  containing  vehicle  or  vemur¬ 
afenib  (0.1  pM  for  A375  cells,  1  pM  for  all  other  cell  lines),  on  dishes  containing 
2.3  X  106  and  6.4  X  106  drug- sensitive  cells,  respectively.  The  same  procedure  was 
followed  for  generation  of  conditioned  media  from  H3 122  (crizotinib,  0.3  or  1  pM) 
or  HCC827  (erlotinib  0.0 1  pM)  lung  adenocarcinoma  cells.  After  72  h,  cells  on  both 
plates  had  reached  equal  confluency  of —80%  and  conditioned  media  was  collected, 
centrifuged  at  1,000  r.p.m.  for  5  min,  filtered,  and  aliquots  were  stored  at  —80  °C 
until  further  use.  Key  proliferation  and  migration  experiments  yielded  the  same 
results  when  performed  with  conditioned  media  in  which  the  same  number  of 
drug-sensitive  cells  (3.2  X  106)  was  plated  initially,  which  resulted  in  higher  cell  con- 
fluency  in  the  vehicle-treated  dish  at  time  of  conditioned  media  collection. 
Proliferation,  survival  and  apoptosis  assays.  Around  1,000-3000  cells  were  pla¬ 
ted  in  a  96-well  plate,  allowed  to  adhere  overnight,  and  then  incubated  with  either 
fresh  or  conditioned  media  containing  vemurafenib  or  additional  drugs  as  indicated. 
After  72  h,  the  number  of  cells  was  determined  using  a  CelltiterGlo  assay  and  the 
caspase  3/7  activity  using  a  CaspaseGlo  assay  (Promega)  according  to  the  manu¬ 
facturer’s  instructions.  Caspase  3/7  activity  was  normalized  to  the  number  of  cells 
present.  All  experiments  with  melanoma  test  cells  and  melanoma  conditioned  media 
were  performed  at  least  three  times,  experiments  with  lung  adenocarcinoma  cell 
lines  were  performed  at  least  twice.  Representative  experiments  are  shown. 
Boyden  chamber  transwell  migration  assay/gap  closure  assay.  Trans  well  migra¬ 
tion  assays  were  performed  as  described  previously  with  minor  modifications31.  In 
brief,  serum-starved  cells  (0.2%  FBS,  overnight)  were  labelled  with  cell  tracker  green 
(Invitrogen)  for  30  min  at  37  °C  and  allowed  to  recover  for  1  h.  Cells  (25,000- 
50,000)  were  then  seeded  onto  membrane  inserts  with  8-pm  pores  and  fluorescence 
blocking  filters  (Falcon).  The  number  of  cells  migrated  through  the  pores  of  the 
membrane  was  scored  after  5-24  h  using  an  Evos  microscope  (AMG).  Gap  closing 
assay  was  performed  according  to  standard  protocols.  In  brief,  cells  were  seeded 
and  grown  until  confluent.  A  tip  was  used  to  generate  a  gap,  cells  were  washed  and 
conditioned  media  was  added.  Images  were  acquired  over  time  to  monitor  for  cap 
closure  in  different  conditions.  All  experiments  were  performed  independently  at 
least  twice.  Representative  experiments  are  shown. 


xCELLigence  migration  assay.  Experiments  were  performed  using  the 
xCELLigence  RTCA  DP  instrument  (Roche  Diagnostics  GmbH)  placed  in  a  humid¬ 
ified  incubator  at  37  °C  with  5%  C02.  Cell  migration  experiments  were  performed 
using  modified  16-well  plates  (CIM-16,  Roche  Diagnostics  GmbH)  according  to 
the  manufacturer’s  instructions.  The  experiment  was  performed  twice.  A  repres¬ 
entative  experiment  is  shown. 

Animal  studies.  All  experiments  using  animals  were  performed  in  accordance  to 
our  protocol  approved  by  MSKCC’s  Institutional  Animal  Care  and  Use  Committee 
(IACUC).  5-7-week-old,  female  NOD-SCID  NCR  (NCI)  or  athymic  NCR-NU- 
NU  (NCI)  mice  were  used  for  animal  experiments  with  human  cell  lines.  Primary 
YUMM1.1  and  YUMM1.7  cell  lines  were  isolated  from  melanomas  developed  in 
mice  ( Tyrr.CreER ;  BrafA ;  Cdkn2a~/~  Ptenlox/lox )  treated  with  4-hydroxytamoxifen 
and  were  subsequently  implanted  in  female  C57BL/6J  (JAX)  mice  aged  between 
5  and  7  weeks.  Tumour  formation,  outgrowth  and  metastasis  were  monitored  by 
BLI  of  TGL-labelled  tumour  cells  as  described  previously22.  In  brief,  anaesthetized 
mice  (150  mg  kg- 1  ketamine,  15  mg  kg- 1  xylazine  or  isoflurane)  were  injected  retro- 
orbitally  with  D-luciferin  ( 1 50  mg  kg- :)  and  imaged  with  an  IVIS  Spectrum  Xenogen 
machine  (Caliper  Life  Sciences).  Bioluminescence  analysis  was  performed  using 
Living  Image  software,  version  4.4.  For  co- implantation  assays,  mice  were  anaes¬ 
thetized  (150  mg  kg-1  ketamine,  15  mg  kg-1  xylazine)  and  1  X  103  TGL-labelled 
resistant  tumour  cells  were  injected  subcutaneously  with  2  X  106  sensitive  tumour 
cells  in  50  pi  growth-factor-reduced  Matrigel/PBS  (1:1)  (BD  Biosciences).  For  the 
control  groups  in  which  the  effects  of  drug  treatment  on  resistant  cells  alone  were 
tested,  2  X  106  resistant  cells  were  injected  in  growth-factor-reduced  Matrigel/ 
PBS.  Two-to-four  sites  on  the  flanks  were  injected  per  mouse.  After  tumours 
reached  a  size  of  50-150  mm3,  the  BLI  signal  of  resistant  cells  was  determined.  To 
compensate  for  minor  growth  differences  of  the  GFP+  resistant  cell  population 
between  mice,  the  mice  were  assigned  to  the  cohorts  so  that  the  overall  BLI  intensity 
(and  consequently  the  cell  number)  was  equal  in  the  treatment  and  control  group. 
Each  group  received  vehicle  or  drug  treatment  as  indicated  (vemurafenib/PLX4032, 
25  mg  kg- 1  twice  daily  for  YUMM 1 . 1  and  YUMM 1 .7  tumours,  and  75  mg  kg- 1  twice 
daily  for  all  other  BRAF  mutant  tumours,  LC-Labs  or  Selleckchem;  100  mg  kg-1 
crizotinib  once  daily,  LC-Labs;  50  mg  kg- 1  erlotinib  once  daily,  LC-Labs;  100  mg  kg- 1 
MK-2206  once  daily,  Chemietek;  50  mg  kg-1  BEZ235  once  daily,  LC-Labs).  Growth 
of  the  resistant  population  in  the  different  groups  was  monitored  by  BLI,  quantified 
and  normalized  to  BLI  signal  at  start  of  treatment.  Tumour  seeding  and  metastasis 
assays  were  performed  as  described  with  minor  modifications20.  In  brief,  sensitive 
tumour  cells  were  injected  subcutaneously  on  two  sites  per  mouse.  Once  tumours 
were  established  (50-150  mm3)  mice  were  treated  with  vehicle  or  vemurafenib 
(75  mg  kg-1  twice  daily)  for  3  days,  and  1  X  105  TGL-labelled  drug- resistant  cells 
were  injected  in  the  left  cardiac  ventricle.  Treatment  was  continued,  and  metastatic 
burden  and  tumour  seeding  were  determined  in  vivo  and  ex  vivo  by  BLI.  Tumour 
volume  was  determined  using  caliper  measurements  and  calculated  using  the  fol¬ 
lowing  formula:  tumour  volume  =  (D  X  d2)/ 2,  in  which  D  and  d  refer  to  the  long 
and  short  tumour  diameter,  respectively.  All  experiments  with  A375  cells  were  in¬ 
dependently  performed  at  least  three  times,  except  animal  experiments  in  Fig.  3, 
which  were  performed  twice.  All  other  animal  experiments  were  independently  per¬ 
formed  at  least  twice.  Representative  experiments  are  shown,  except  where  noted 
and  where  instead  the  average  of  three  experiments  is  presented. 

Gene  expression  analysis.  Whole  RNA  was  isolated  from  cells  using  RNAeasy  Mini 
Kit  (QIAGEN).  The  Transcriptor  First  Strand  cDNA  synthesis  kit  (Roche)  was  used 
to  generate  cDNA.  Differential  RNA  levels  were  assessed  using  Taqman  gene 
expression  assays  (Life  technologies).  Assays  used  for  human  genes  are: 
Hs041 87685,  Hs00365742,  Hs00605382,  Hs00601975,  Hs01099999,  Hs00959010, 
HsO  1029057,  Hs00234244,  Hs00905117,  Hs00180842,  Hs00989373,  Hs00234140, 
Hs00195591,  Hs00207691,  Hs99999141,  HsOl  117294,  Mm00607939,  Mm99999915 
and  Mm04207958.  Relative  gene  expression  was  normalized  to  internal  control  genes: 
B2M  (Hs99999907_ml),  GAPDH  (Hs99999905_ml)  zndACTB  (Mm00607939_ 
si).  Quantitative  PCR  reactions  were  performed  on  a  VIIA7  Real-Time  PCR  system 
and  analysed  using  VII A7  software  (Life  Technologies).  All  data  points  represent  at 
least  four  technical  replicates  and  experiments  were  performed  independently  three 
times.  A  representative  experiment  is  shown. 

Cancer-cell-specific  TRAP  and  sequencing.  To  investigate  the  gene  expression 
changes  specifically  of  drug-sensitive  tumours  during  vemurafenib  treatment,  or  gene 
expression  changes  of  resistant  cells  exposed  to  a  regressing  tumour  microenvi¬ 
ronment,  A375  and  A375R  cells,  respectively,  were  modified  to  express  eGFP-RPLIOa. 
Tumours  derived  from  implanted  A375-eGFP-RPL10a  and  A375R-eGFP-RPL10a 
cells  were  homogenized  and  processed  with  the  TRAP  protocol  as  previously 
described28,32,33  with  the  following  modifications:  fresh  tumour  was  homogenized 
with  a  Model  PRO  200  homogenizer  at  speed  5  for  four  cycles  of  15  s,  RNasin  Plus 
RNase  inhibitor  (Promega,  N2615)  was  used  as  RNase  inhibitor,  and  anti-eGFP 
antibody  coated  sepharose  beads  (GE  Healthcare)  were  used  for  immunopreci- 
pitation.  Polysome-associated  RNA  was  purified  with  RNAqueous  micro  kit  (Life 
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Technologies,  AM1931).  Ribogreen  and  the  Agilent  BioAnalyzer  technologies  were 
used  to  quantify  and  control  the  quality  of  RNA;  500  ng  RNA  (RNA  integrity  num¬ 
ber  (RIN)  >  8.5)  from  each  sample  was  used  for  library  construction  with  TruSeq 
RNA  Sample  Prep  Kit  v2  (Illumina)  according  to  the  manufacturer’s  instructions. 
The  samples  were  barcoded  and  run  on  a  Hiseq  2000  platform  in  a  50-base-pair 
(bp)/50-bp  or  75-bp/75-bp  paired-end  run,  using  the  TruSeq  SBS  Kit  v3  (Illumina). 
An  average  of  40  million  paired  reads  was  generated  per  sample. 

RNA-seq  analysis.  For  drug- sensitive  A375,  C0I08OO,  UACC63  and  H3122  cells, 
in  vitro ,  raw  paired-end  sequencing  reads  were  mapped  to  the  human  genome  (build 
hgl9)  with  STAR2.3.0e  (ref.  34)  using  standard  options.  Uniquely  mapped  reads 
were  counted  for  each  gene  using  HTSeq  vO.5.4  (ref.  35)  with  default  settings.  Read 
counts  of  each  sample  were  normalized  by  library  size  using  the  cDESeq’35  package 
of  Bioconductor.  Differential  gene  expression  analysis  between  any  two  conditions 
was  performed  based  on  a  model  using  the  negative  binomial  distribution35.  Genes 
with  false  discovery  rate  (FDR)  <  0.05,  fold  change  larger  than  1.5  or  smaller  than 
0. 667-fold,  and  average  read  counts  larger  than  10  were  treated  as  differentially 
expressed  genes.  RNA-seq  data  from  in  vivo  xenograft  TRAP  samples  were  pro¬ 
cessed  with  the  following  modifications  to  avoid  potential  mRNA  contamination 
from  host  mouse  tissue:  raw  sequencing  reads  were  mapped  to  a  hybrid  genome 
consisting  indexes  of  both  human  (build  hgl9)  and  mouse  (build  mm9)  genomes. 
Only  reads  that  uniquely  mapped  to  human  genome  indexes  were  preserved  and 
counted  using  HTSeq  vO.5.4  (ref.  35). 

Bioinformatics  analysis.  Heatmap  visualization  of  data  matrices  was  performed 
using  the  ‘gplots’  package  of  R.  Principle  component  analysis  of  RNA-seq  results  was 
performed  with  the  variance  stabilizing  transformation  methods  in  cDESeq’  package 
of  Bioconductor  and  the  first  two  principal  components  were  plotted.  Volcano  plots 
were  derived  from  cDESeq’-based  differential  gene  expression  analysis.  Differentially 
expressed  genes  with  transcription  factor  activity  (G0:00037000)  at  6  h  ofvemur- 
afenib  treatment  and  gene  products  located  in  the  extracellular  region  (GO: 
00005576)  at  48  h  of  vemurafenib  treatment  were  identified  using  the  Database  for 
Annotation,  Visualization  and  Integrated  Discovery  (DAVID)36  v6.7  (http://david. 
abcc.ncifcrf.gov/)  and  enriched  GO  terms  were  visualized  using  RE  VIGO37  (http:// 
revigo.irb.hr).  Enriched  transcriptional  regulators  for  the  list  of  differentially  ex¬ 
pressed  gene  products  in  the  extracellular  region  were  predicted  with  DAVID  v6.7 
and  this  list  compared  to  the  gene  expression  levels  of  transcription  factors  after 
6  h  of  vemurafenib  treatment  in  A375  cells.  Upstream  regulators,  functions  assoc¬ 
iated  with  the  gene  expression  profile  and  potential  drug  vulnerabilities  were  deter¬ 
mined  by  interpretative  phenomenological  analysis  (IPA)  analysis  on  differentially 
expressed  genes  from  A375R-eGFP-RPL10a  cells  in  different  tumour  microenvir¬ 
onments  in  vivo. 

Immunoblotting.  RIPA  buffer  (Cell  Signaling)  was  used  for  cell  lysis,  according  to 
the  manufacturer’s  instructions,  and  the  protein  concentrations  were  determined 
by  BCA  Protein  Assay  kit  (Pierce).  Proteins  were  separated  by  SDS-PAGE  using 
Bis-Tris  4-12%  gradient  polyacrylamide  gels  in  the  MOPS  buffer  system  (Invitro- 
gen)  and  transferred  to  nitrocellulose  membranes  (Bio Rad)  according  to  standard 
protocols.  Membranes  were  immunoblotted  with  antibodies  against  pERKT202/Y204 
(4370),  tERK  (4696),  pAKTS473  (4060),  pAKTX308  (40  56),  tAKT  (2920),  EGFR 
(4267),  MET  (8198),  PDGFRb  (3169),  pFRAl  (3880),  caspase3  (9662),  pPRAS40T246 
(13175),  p70S6KT389  (92  05),  pFAKY397  (32  83),  pPKCbetaIIS660  (93  71),  pNFKBS536 
(3033),  pP-CateninS33/37/T41  (9561),  pSTAT-3Y705  (9  1  45),  pSTAT-5CUC5  (9359), 
pGSK3a/pS21/9  (9327),  pCREBsl33/pATF-l  (9196)  (Cell  Signaling,  1:1,000),  FRA1 
(sc-605,  Santa  Cruz  Biotechnology,  1:200)  and  tubulin  (T6074,  Sigma-Aldrich, 
1 :5,000)  in  Odyssey  blocking  buffer  (LI-COR) .  After  primary  antibody  incubation, 
membranes  were  probed  with  IRDye  800CW  donkey- anti- mouse  IgG  (LI-COR) 
or  IRDye  680RD  goat-anti-rabbit  IgG  (LI-COR)  secondary  antibody  (1:20,000) 
and  imaged  using  the  LI-COR  Odyssey  system.  All  immunoblots  were  performed 
independently  at  least  twice.  Tubulin  served  as  a  loading  control. 

Plasmids,  recombinant  protein  and  ELISA.  Identifiers  for  shRNAs  used  in  this 
study  are:  V3LHS-644610  (shFRAl-1),  V3LHS-644611  (shFRAl-2),  V3LHS-320021 
(shIGFBP3-l)  and  V2LHS- 1 1 1629  (shIGFBP3-2)  (Dharmacon,  GELifesciences). 
IGFBP3  ELISA  (Raybiotech)  was  performed  according  to  the  manufacturer’s  in¬ 
structions  with  50  pg  tumour  lysate  and  conditioned  media  was  diluted  1 :5.  Recom¬ 
binant  proteins  were  used  at  the  following  conditions:  10  ng  ml-1  IGF1  (Invitrogen), 
lOngml-1  EGF  (Invitrogen),  lOngml-1  PDGFD  (R&D  Systems),  2pgml-1 
IGFBP3  (Prospec)  for  15  min,  or  5  pg  ml-1  ANGPTL7  (R&D  Systems)  for  30  min. 
Patient  samples.  Melanoma  tissues  were  obtained  from  clinical  trial  patients  or 
patients  under  standard  clinical  management  with  approval  of  the  UCLA  Institu¬ 
tional  Review  Board.  Patient-informed  consent  was  obtained  for  the  research  per¬ 
formed  in  this  study. 

Immunofluorescence.  Tissues  for  BrdU-immunofluorescence  staining  were  ob¬ 
tained  after  overnight  fixation  with  4%  paraformaldehyde  (PFA)  at  4  °C,  embed¬ 
ded  in  OCT  compound  (VWR)  and  stored  at  —80  °C.  10-pm  thick  cryosections 
on  glass  slides  were  used  for  immunofluorescence  staining  according  to  standard 


protocols.  Tissue  for  all  other  immunofluorescence  experiments  from  xenograft 
tumours  was  obtained  after  fixation  with  4%  PFA  at  4  °C  and  a  series  of  dehydra¬ 
tion  steps  from  15%  to  30%  sucrose,  as  described  previously38.  In  brief,  tumours 
were  sliced  using  a  sliding  microtome  (Fisher).  Tumour  slices  (80  pm)  were  blocked 
floating  in  10%  NGS,  2%  BSA,  0.25%  Triton  in  PBS  for  2  h  at  room  temperature. 
Primary  antibodies  were  incubated  overnight  at  4  °C  in  the  blocking  solution  and 
the  next  day  for  30  min  at  room  temperature.  After  washes  in  PBS-Triton  0.25%, 
secondary  antibodies  were  added  in  the  blocking  solution  and  incubated  for  2  h. 
After  extensive  washing  in  PBS-Triton  0.25%,  nuclei  were  stained  with  Bis- 
Benzamide  for  5  min  at  room  temperature,  tumour  slices  were  washed  and  trans¬ 
ferred  to  glass  slides.  Slices  were  mounted  with  ProLong  Gold  anti  fade  reagent 
(Invitrogen).  Primary  antibodies:  GFP  (GFP-1020,  Aves  Labs,  1:1,000),  collagen  IV 
(AP756,  Millipore,  1:500),  BrdU  (ab6326,  Abeam,  1:250),  FRA1  (sc605,  Santa  Cruz, 
1:200).  Secondary  antibodies:  Alexa-Fluor-488  anti-chicken,  Alexa-Fluor-555  anti¬ 
rabbit,  Alexa-Fluor-555  anti-rat  (Invitrogen).  Stained  sections  were  visualized  using 
a  Carl  Zeiss  Axioimager  Z1  microscope  or  with  a  Leica  SP5  upright  confocal  micro¬ 
scope  using  X 10  or  X  20  objectives.  Images  were  analysed  with  Image J,  and  Meta- 
morph  software. 

Flow  cytometry.  Flow  cytometry  was  performed  as  described  previously22,  with 
minor  modifications.  In  brief,  whole  tumours  were  dissected,  cut  into  smaller  sec¬ 
tions  and  dissociated  for  1-3  h  with  0.5%  collagenase  type  III  (Worthington  Bio¬ 
chemical)  and  1%  dispase  II  (Roche)  in  PBS.  Resulting  single  cells  suspensions  were 
washed  with  PBS  supplemented  with  2%  FBS  and  filtered  through  a  70-pm  nylon 
mesh.  The  resulting  single  cell  suspension  was  incubated  for  10  min  at  4  °C  with 
anti-mouse  Fc-block  CD  16/32  antibody  (2.4G2  BD)  in  PBS  supplemented  with 
1%  BSA.  Cells  were  subsequently  washed  with  PBS/BSA  and  stained  with  control 
antibodies  or  antibodies  to  detect  immune  cells  diluted  in  PBS  supplemented  with 
0.5%  BSA  and  2  mM  EDTA.  The  following  antibodies  against  mouse  antigens  were 
used:  CD45-PE-Cy7  (clone  30-F11,  BD  Pharmingen,  1:200),  CDllb-APC  (clone: 
Ml/70,  BD  Pharmingen,  1:100),  Grl-PE  (MACS,  1:10),  CD31-APC  (clone:  390, 
eBioscience,  1:100),  F4/80-PE  (clone:  BM8,  eBioscience,  1:50).  To  determine  the 
level  of  EdU  incorporation  in  A375R  cells  within  vehicle-  or  vemurafenib-treated 
A375/A375R  tumours,  EdU  (50  mg  kg-1,  Life  Technologies)  was  injected  intraper- 
itoneally,  after  2  h  tumours  were  collected,  single-cell  suspensions  generated  as 
described  above  and  further  processed  according  to  the  manufacturer’s  protocol 
(Click- iT  Plus  EdU  Alexa  Fluor  647  Flow  Cytometry  Assay  Kit,  Life  technologies). 
Data  were  acquired  using  a  FACS  Calibur  (BD  Biosciences).  All  experiments  were 
performed  independently  at  least  two  times.  Representative  experiments  are  shown. 
Antibody  arrays.  Cytokines  and  cytokine  receptors  of  murine  stromal  and  immune 
cells,  in  A375  tumours  treated  with  vehicle  or  vemurafenib  for  5  days,  were  mea¬ 
sured  using  the  Mouse  Cytokine  Array  G2000  (RayBio,  AAH-CYT-G2000-8,  de¬ 
tecting  174  proteins),  according  to  the  recommended  protocols.  In  brief,  tumours 
were  homogenized  with  a  Mini  Immersion  Blender  (Pro  Scientific)  in  Raybio  Lysis 
buffer  with  protease  inhibitors.  Lysates  were  centrifuged  for  5  min  at  10,000g,  the 
supernatant  was  collected  and  protein  concentration  was  measured  using  the  BCA 
Assay  Kit  (Pierce).  Protein  (150  pg)  was  hybridized  on  the  antibody  arrays  over¬ 
night  at  4  °C.  IRDye-labelled  streptavidin  (LI-COR)  at  a  dilution  of  1 :5,000  was  used 
for  the  detection,  slides  were  scanned  using  a  Odyssey  CLx  scanner  (LI-COR)  and 
analysed  using  Image  Studio  2.0  software.  The  results  were  then  normalized  using 
internal  controls,  and  the  relative  protein  levels  determined  across  four  biological 
replicates. 

Senescence  P-galactosidase  staining.  A375  cells  were  grown  in  low-serum  media 
and  treated  with  vehicle  or  vemurafenib  (0.1  pM)  for  3  or  8  days,  P-galactosidase 
staining  was  performed  according  to  the  manufacturer’s  instructions  (Cell  Signaling). 
All  experiments  were  performed  independently  three  times.  Representative  ex¬ 
periments  are  shown. 

Statistical  analysis.  Data  are  generally  expressed  as  mean  ±  s.e.m.,  or  in  box  plots 
in  which  the  centre  line  is  the  median,  and  whiskers  are  minimum  to  maximum 
values.  Group  sizes  were  determined  based  on  the  results  of  preliminary  experiments 
and  no  statistical  method  was  used  to  predetermine  sample  size.  Group  allocation 
and  outcome  assessment  were  not  performed  in  a  blinded  manner.  All  samples  that 
met  proper  experimental  conditions  were  included  in  the  analysis.  Statistical  sig¬ 
nificance  was  determined  using  a  two-tailed  Mann- Whitney  U  test  or  Student’s 
f-test  using  Prism  6  software  (GraphPad  Software),  or  using  a  hypergeometric  var¬ 
iability  test  (http://www.geneprof.org).  Significance  was  set  at  P  <  0.05. 
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Extended  Data  Figure  1  |  Targeted  therapy  or  oncogene  knockdown  leads  to 
regression  of  sensitive  melanoma  and  lung  adenocarcinoma  tumours  but 
accelerates  the  proliferation  and  seeding  of  residual  drug-resistant  cells 

in  vivo,  a,  FACS  analysis  of  sensitive  A375  and  vemurafenib- resistant  A375R 
cells  expressing  TGL,  at  tumour  implantation  and  after  2  weeks  at  start  of 
therapy  (n  —  8  tumours).  Plots  depict  representative  images,  b,  Tumour 
volume  of  A375  cells  treated  with  vehicle  or  vemurafenib  over  time  (vehicle, 
n  =  8;  vemurafenib,  n  —  12  tumours),  c,  Representative  sections  of  A375/ 
A375r-TGL  tumours  at  0,  1,3  and  6  days  of  vemurafenib  treatment  analysed 
with  immunofluorescence  against  GFP.  Arrowheads  indicate  emerging  clusters 
of  GFP+  resistant  cells.  Scale  bars,  2  mm.  d,  Quantification  of  BrdU 
incorporation  into  vemurafenib -resistant  A375R-TGL  cells  in  A375/A375R 
tumours  treated  with  vehicle  or  vemurafenib  for  6  days  (vehicle,  n  —  13  FOV  of 
3  tumours;  vemurafenib  n  =  18  FOV  of  4  tumours).  Original  magnification, 
X20.  e,  Fold  change  of  photon  flux  of  TGL- expressing  A375R  cells  in  A375 
tumours  or  A375R  tumours  alone  treated  with  vehicle  or  dabrafenib  for  8  days 
(A375/A375r:  vehicle,  n  =  15;  dabrafenib,  n  =  14;  A375R:  vehicle,  n  =  8; 
dabrafenib,  n  —  7  tumours),  f,  Tumour  volume  of  doxycycline-inducible  BRAF 
knockdown  A375-i-shBRAF-derived  xenograft  tumours  (in  which  ‘A375-i’ 
denotes  expression  of  doxycycline-inducible  hairpin)  treated  with  vehicle  or 
doxycycline  over  time  (vehicle,  n  —  5;  doxycycline,  n  —  4  tumours),  g,  Photon 
flux  of  TGL- expressing  A375R  cells  mixed  in  A375-i-shBRAF  tumours  treated 
with  vehicle  or  doxycycline  (vehicle,  n  —  10;  doxycycline,  n—  11  tumours), 
h,  Fold  change  of  photon  flux  of  TGL-expressing  vemurafenib-resistant 
M249R4  tumours  treated  with  vehicle  or  vemurafenib  (n  =  16  tumours). 


i-k,  Co-implantation  assay  of  tumours  treated  with  vehicle  or  corresponding 
targeted  therapy  with  BLI  quantification  after  5-8  days,  i,  Fold  change  of 
photon  flux  of  TGL-expressing  vemurafenib-resistant  YUMM1.7R  cells  mixed 
in  unlabelled,  vemurafenib-sensitive  YUMM1.7  tumours  or  YUMM1.7R 
tumours  alone  (YUMM1.7/YUMM1.7R:  n  =  24  ;  YUMM1.7R:  n  =  20 
tumours),  j,  Fold  change  of  photon  flux  of  TGL-expressing,  intrinsically 
vemurafenib-resistant  B16  cells  mixed  in  vemurafenib-sensitive  YUMM1.1 
tumours  or  B16  tumours  alone  (YUMM1.1/B16:  vehicle,  n  =  12;  vemurafenib, 
n  —  16;  B16:  n  —  20  tumours),  k,  A375R  mixed  in  crizotinib- sensitive  H3122 
cells  or  A375r  tumours  alone  (H3122/A375R:  vehicle,  n  =  14;  crizotinib, 
n  =  13;  A375r:  n  —  12  tumours).  1,  Photon  flux  of  tumours  established  from 
intrinsically  resistant  drug-resistant  cells  alone,  treated  with  vehicle,  crizotinib 
or  erlotinib  (crizotinib-resistant  PC9,  H2030-  or  erlotinib-resistant  A375R) 

( n  (from  left  to  right)  =  12,  12,  7,  12,  16  and  16  tumours,  respectively), 
m,  Summary  of  the  model  systems  and  conditions  used  in  vivo,  n,  Left, 
representative  immunofluorescence  images  of  vemurafenib-treated,  sensitive 
tumours  7  h  or  5  days  after  intracardiac  injection  with  A375R-TGL  cells; 
sections  stained  for  GFP  (A375R,  green),  collagen  type  IV  (blood  vessels,  red), 
and  DAPI  (nuclei,  blue).  Right,  quantification  of  A375R  single  cells  and  cell 
clusters  (>2  cells)  infiltrating  an  A375  tumour  treated  with  vehicle  or 
vemurafenib  after  intracardiac  injection  of  A375R  cells  (GFP+  cells  were  scored 
in  at  least  10  whole  sections  of  at  least  4  tumours).  Original  magnifications, 
X20.  Data  in  b,  e-1  and  n  are  mean  and  s.e.m.,  in  f,  centre  line  is  median, 
whiskers  are  minimum  to  maximum.  P  values  calculated  by  a  two-tailed 
Mann-Whitney  U  test. 
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Extended  Data  Figure  2  |  The  secretome  of  vemurafenib-treated  melanoma 
and  crizotinib-  or  erlotinib-treated  lung  adenocarcinoma  cells  stimulates 
the  proliferation  and  migration  of  drug-resistant  cells  in  vitro  and  occurs 
before  apoptosis  and  senescence,  a,  Quantification  of  the  co-culture  assay, 
depicted  in  Fig.  2a,  7  days  after  addition  of  resistant  A375R-TGL  cells  (n  =  4 
biological  replicates).  P  values  calculated  using  a  Student’s  t-test.  b,  c,  Drug- 
sensitive  cells  were  pre-treated  with  vehicle  or  drug  (crizotinib  or  erlotinib) 
for  48  h  before  5  X  102  TGL-expressing,  drug-resistant  cells  were  added. 
Growth  was  monitored  by  BLI  and  quantified  7  days  after  addition  of  the 
resistant  cell  population  (n  —  8  biological  replicates),  P  values  calculated  using  a 
Student’s  t-test.  b,  Quantification  and  representative  images  of  TGL-expressing 
H2030  cells  alone  or  co-cultured  with  crizotinib- sensitive  H3122  cells  and 
treated  with  vehicle  or  crizotinib  c,  Quantification  and  representative  images  of 
TGL-expressing  A375R  cells  alone  or  co-cultured  with  erlotinib -sensitive 
HCC827  cells  and  treated  with  vehicle  or  erlotinib.  d,  Relative  number  of 
vemurafenib- resistant  LOXR  cells  after  3  days  in  the  presence  of  conditioned 
media  derived  from  A375  and  UACC62  cells  (n  =  3  biological  replicates), 
e,  Representative  immunofluorescence  for  Ki67  in  drug-resistant  YUMM1.7R 
cells  cultured  in  conditioned  media  from  YUMM1.7  cells.  Original 
magnification,  X20.  f,  Relative  number  of  vemurafenib- resistant  melanoma 
cells  with  different,  clinically  relevant  resistance  mechanisms  after  3  days  in  the 
presence  of  conditioned  media  derived  from  A375  cells.  SKMEL239-3 
expressing  the  p61  BRAFV600E  splice  variant,  A375  expressing  NRASQ61K  or 


the  constitutively  active  MEK  variant  MEK-DD  (n  —  5  biological  replicates), 
g,  Relative  cell  number  of  intrinsically  vemurafenib -resistant  lung 
adenocarcioma  cells  (H2030,  PC9)  or  crizotinib-  and  erlotinib -resistant 
melanoma  cells  (A375R)  after  3  days  cultured  in  the  presence  of  conditioned 
media  from  vemurafenib-treated  melanoma  or  crizotinib-  and  erlotinib- 
treated  lung  adenocarcinoma  (n  =  6  in  all,  except  for  A375R  with  HCC827- 
CM,  n  —  4  biological  replicates),  h,  Representative  image  of  A375R  cells 
migrated  towards  A375-derived  CM-vehicle  or  CM -vemurafenib.  Original 
magnification,  X10.  i,  Relative  migration  towards  conditioned  media 
from  different  sources  and  different  resistant  test  cells  as  indicated  (n  —  10 
FOV).  **P<0.01,  ****p<  0.0001,  two-tailed  Mann- Whitney  L/test. 
j,  Representative  graph  and  quantification  of  real-time  migration  of  A375R  cells 
in  the  presence  of  conditioned  media  derived  from  A375  cells  as  measured 
by  the  xCELLigence  system  {n  —  4  biological  replicates).  Rvalue  calculated 
using  two-tailed  Mann- Whitney  U  test,  k,  Monolayer  gap  closing  assay  of 
A375r  cells  in  the  presence  of  conditioned  media  derived  from  A375  cells  with 
representative  light  microscopy  images  and  quantification  of  gap  closure  over 
time.  1,  Immunoblotting  for  cleaved  caspase-3  and  phosphorylated  ERK 
protein  levels  in  vemurafenib -sensitive  melanoma  cell  lines  after  72  h  of 
vemurafenib  treatment,  m,  (3-galactosidase  staining  of  A375  cells  treated  with 
vemurafenib  for  72  h  or  8  days.  Original  magnification,  X 13.  Data  are 
presented  mean  and  s.e.m. 
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Extended  Data  Figure  3  |  The  therapy-induced  secretome  of  sensitive  cells 
overlaps  significantly  in  melanoma  and  lung  adenocarcinoma  cells  and 
appears  after  gene  expression  changes  enriched  for  transcriptional 
regulators,  a,  b,  GO  analysis  (http://revigo.irb.hr)  of  gene  expression  changes 
after  6  h  of  vemurafenib  treatment  of  A375  cells  with  spatial  representation  of 
enriched  GO  terms  (a)  and  the  molecular  functions  significantly  affected  (b). 
c,  Heat  map  representing  the  expression  levels  of  commonly  up-  and 
downregulated  genes  in  vemurafenib-treated  A375-derived  xenograft  tumours 
(5  days)  and  A375  cells  in  vitro  (48  h).  d,  Principal  component  analysis  of 
vemurafenib- sensitive  C0I08OO  and  UACC62  melanoma  cells  and  crizotinib- 
sensitive  H3122  lung  adenocarcinoma  cells  treated  in  vitro  with  vehicle  or 
vemurafenib  or  crizotinib  for  48  h.  e,  Venn  diagram  indicating  the  overlap  of 
genes  in  the  extracellular  region  (G0:0005576)  upregulated  after  48  h  of 


vemurafenib  treatment  in  A375,  C0I08OO  and  UACC62  melanoma  cell  lines, 
f,  Venn  diagram  indicating  the  overlap  of  genes  in  the  extracellular  region 
(G0:0005576)  upregulated  after  48  h  of  vemurafenib  treatment  in  at  least  2  out 
of  3  melanoma  models  and  after  48  h  of  crizotinib  treatment  in  the  H3122 
lung  adenocarcinoma  cell  line,  g,  Venn  diagram  indicating  the  overlap  of  genes 
in  the  extracellular  region  (G0:0005576)  downregulated  after  48  h  of 
vemurafenib  treatment  in  A375,  C0I08OO  and  UACC62  melanoma  cell  lines, 
h,  Venn  diagram  indicating  the  overlap  of  genes  in  the  extracellular  region 
(G0:0005576)  downregulated  after  48  h  of  vemurafenib  treatment  in  at  least 
2  out  of  3  melanoma  models  and  after  48  h  of  crizotinib  treatment  in  the  H3 122 
lung  adenocarcinoma  cell  line.  P  values  calculated  using  a  hypergeometric 
probability  test. 
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Extended  Data  Figure  4  |  Vemurafenib  treatment  induces  widespread 
changes  in  the  intra-tumour  immune  cell  composition  and  stromal  cytokine 
composition  in  tumours  regressing  during  targeted  therapy,  a,  b,  FACS 
analysis  of  murine  immune  cell  populations  in  A375-derived  xenograft 
tumours  treated  with  vehicle  or  vemurafenib  for  5  days,  a,  b,  Representative 
image  (a)  and  quantification  (b)  of  intra-tumour  composition  of  indicated 


immune  cell  populations  (vehicle,  n  —  4;  vemurafenib,  n  —  6  tumours), 
c,  d,  Cytokine  array  of  murine  stroma- derived  cytokines  within  A375-derived 
xenograft  tumours  treated  with  vehicle  or  vemurafenib  for  5  days. 
Representative  image  (c)  and  quantification  (d)  of  down-  and  upregulated 
cytokines  {n  —  4  tumours).  P  values  calculated  by  a  two-tailed  Mann- Whitney 
U  test.  Data  are  mean  and  s.e.m. 
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Extended  Data  Figure  5  |  Targeted  therapy  induces  downregulation  of 
FRA1  in  drug-sensitive  tumour  cells,  a,  List  of  transcription  factors  predicted 
to  regulate  the  vemurafenib-induced  reactive  secretome  in  A375  cells,  and  a 
heat  map  of  the  corresponding  transcription  factor  gene  expression  levels 
in  these  cells.  Red  represents  high,  yellow  medium  and  blue  low  relative 
expression  on  the  colour  scale,  b,  Immunoblotting  of  phosphorylated  and  total 
FRA1  protein  levels  in  A375  and  YUMM1.7  melanoma  cell  lines  treated 
with  vemurafenib  for  24  h.  c,  Relative  mRNA  levels  of  FRA1  in  H3122  cells 


treated  with  crizotinib  (500  nM)  and  HCC827  treated  with  erlotinib  (10  nM)  at 
different  time  points  {n  —  4  technical  replicates),  d,  Relative  mRNA  levels 
of  FRA1  in  A375R  cells  treated  with  vemurafenib  at  different  time  points  {n  —  4 
technical  replicates),  e,  Immunofluorescence  staining  of  FRA1  (red)  and 
DAPI  (blue)  in  biopsies  from  melanoma  patients  before  and  after  vemurafenib 
treatment  (clinical  information  can  be  found  in  Extended  Data  Table  1). 
Original  magnification,  X20. 
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Extended  Data  Figure  6  |  The  secretome  of  melanoma  cells  with  FRA1 
knockdown  stimulates  proliferation  and  migration  of  A375R  cells  in  vitro 
and  in  vivo,  a,  Immunoblotting  of  phosphorylated  and  total  FRA1  protein 
levels  in  A375  cells  transduced  with  control  shRNA,  with  or  without  additional 
vemurafenib  treatment,  or  shRNAs  targeting  FRA1.  b,  Photon  flux  and 
representative  BLI  images  of  TGL- expressing  A375R  cells  co-cultured  with 
A375  cells  expressing  control  shRNA  (with  or  without  vemurafenib  treatment) 
or  FRA  1  -targeting  shRNAs  after  7  days  (n  —  9  biological  replicates),  c,  Relative 
number  of  A375R  cells  after  3  days  in  the  presence  of  conditioned  media  derived 
from  A375  cells  transduced  with  control  shRNA,  with  or  without  additional 
vemurafenib  treatment,  or  FRA1  shRNAs  (n  —  3  biological  replicates). 


d,  Migration  of  A375R  cells  towards  conditioned  media  derived  from  A375  cells 
transduced  with  control  shRNA  (with  or  without  vemurafenib  treatment)  or 
FRA1  shRNAs  using  a  Boyden  chamber  assay  (shCtrl,  n  =  15;  all  other  groups 
n  —  10  FOV)  e,  Relative  mRNA  levels  of  selected  secreted  factors  and 
transcription  factors  of  A375  cells  expressing  control  shRNA  or  an  shRNA 
targeting  FRA1  (shFRAl-1),  treated  with  vehicle  or  vemurafenib  (24  h). 
f,  Bioluminescent  signal  of  A375R-TGL  cells  8  days  after  subcutaneous  co¬ 
implantation  with  UACC62  cells  expressing  a  control  or  an  shRNA  for  FRA1 
(shCtrl,  n  =  12;  shFRAl,  n  =  20  tumours).  Data  are  mean  and  s.e.m.  *P  <  0.05, 
**P  <  0.01,  ***P  <  0.001,  ****P  <  0.0001,  Student’s  t-test. 
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Extended  Data  Figure  7  |  The  TIS  includes  upregulated  positive  regulators 
and  a  loss  of  negative  regulators  of  the  PI(3)K/AKT/mTOR  pathway, 
which  is  activated  in  sensitive  and  resistant  cells  in  vitro  and  in  vivo. 

a,  b,  Enriched  biological  processes  (a)  and  inferred  drug  vulnerabilities  (b)  as 
determined  by  Ingenuity  pathway  analysis  of  gene  expression  data  from 
vemurafenib-resistant  A375R  cells  responding  to  signals  from  the  reactive 
tumour  microenvironment  of  a  tumour  regressing  during  targeted  therapy 
in  vivo  (for  experimental  set-up  see  Fig.  la  and  Methods),  c,  Left, 
immunoblotting  of  phosphorylated  AKTS473  and  phosphorylated  ERK  protein 
levels  in  A3 75  cells  treated  with  vehicle  or  vemurafenib  at  different  time  points 
during  the  generation  of  conditioned  media.  Right,  immunoblotting  of 
phosphorylated  AKTS473  and  phosphorylated  ERK  protein  levels  in  A375 
cells  after  short-term  exposure  to  conditioned  media  derived  from  A375  cells 
treated  with  vehicle  or  vemurafenib.  d,  Immunoblotting  of  phosphorylated 
AKTS473  and  phosphorylated  FRA1  protein  levels  in  A375-derived  xenograft 
tumours  treated  with  vehicle  or  vemurafenib  for  5  days.  Normalized 
quantification  of  phospho-AKTS473/tubulin  in  the  bottom  panel, 
e,  Immunoblotting  of  a  range  of  pathway  nodes  in  A375R  cells  treated  with 


CM-vehicle  or  CM -vemurafenib,  derived  from  A375  cells,  for  15,  30,  60  or 
120  min.  f,  Cancer  cell-derived  IGFBP3  levels  (left)  and  murine  stromal  IGF1 
levels  (right)  in  A375-derived  xenograft  tumours  treated  with  vehicle  or 
vemurafenib  for  5  days  as  determined  by  ELISA  {n  —  4  tumours)  g,  Cancer¬ 
cell-derived  IGFBP3  levels  in  conditioned  media  from  indicated  melanoma  cell 
lines  treated  with  vehicle  or  vemurafenib  as  determined  by  ELISA  (n  —  3 
technical  replicates  of  conditioned  media  derived  from  at  least  two  biological 
replicates),  h,  IGFBP3  levels  in  conditioned  media  derived  from  A375  cells 
expressing  control  shRNA  or  shRNAs  targeting  IGFBP3  (shIGFBP3-l  and  -2) 
as  determined  by  ELISA  (n  =  3  technical  replicates),  i,  Immunoblotting  of 
phosphorylated  AKTS473  in  A375R  cells  after  incubation  with  conditioned 
media  of  A375  cells  expressing  control  shRNA  or  shRNAs  targeting  IGFBP3. 
j,  Phosphorylation  status  of  AKTS473  in  A375R  cells  after  incubation  for  15  min 
with  conditioned  media,  IGF1  and  IGFBP3  as  indicated,  k,  Bioluminescent 
signal  of  A375R-TGL  cells  10  days  after  co-implantation  with  A375  cells 
expressing  a  control  shRNA  or  an  shRNA  targeting  IGFBP3  (shIGFBP3-l) 

{n  —  10  tumours).  P  values  calculated  by  a  two-tailed  Mann- Whitney  U  test. 
Data  are  mean  and  s.e.m. 
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Extended  Data  Figure  8  |  Dual  inhibition  of  RAF  and  the  AKT/mTOR 
pathway  blunts  the  effects  of  the  regressing  tumour  environment  on  the 
resistant  cell  population,  a,  Relative  photon  flux  and  representative  BLI 
images  of  GFP/luciferase  expressing  A375R  cells  co-cultured  with  A375  cells 
and  treated  with  vehicle,  vemurafenib  or  the  combination  of  vemurafenib  and 
either  MK2206  (AKTi,  2  pM)  or  BEZ235  (PI(3)K/mTORi,  300  nM)  for  7  days 
( n  —  2-3  biological  replicates),  b,  Relative  number  of  A375R  cells  after  3  days 
in  the  presence  of  CM-vehicle  or  CM -vemurafenib  with  additional  BEZ235 
(300  nM)  (n  =  3  biological  replicates),  c,  Mice  bearing  tumours  consisting  of 
A375/A375R  cells  or  A375R  cells  alone  were  treated  with  drugs  as  indicated. 


Bioluminescent  signal  of  TGL-expressing  A375R  cells  was  determined  on  day  5 
of  treatment  ( n  =  16,  16,  12,  12,  12  and  16  tumours,  respectively),  d,  Mice 
bearing  tumours  consisting  of  unlabelled  A375  cells  were  pre-treated  for  3  days 
with  drugs  as  indicated  and  1  X  105  TGL-expressing  A375R  cells  were 
inoculated  in  the  arterial  circulation.  Drug  treatment  was  continued  and 
seeding  of  resistant  cells  to  the  primary  tumour  was  quantified  by  BLI. 
Representative  BLI  images  on  the  right  (vehicle,  n  —  4;  vemurafenib  n  =  10, 
vemurafenib  plus  BEZ235,  n  —  10  tumours).  P  values  calculated  by  a  two-tailed 
Mann- Whitney  test.  Data  are  mean  and  s.e.m. 
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Extended  Data  Figure  9  |  Characterization  of  cell  lines  in  response  to 
targeted  therapy,  a-h,  Relative  survival  of  human  melanoma  cell  lines 
(A375,  C0I08OO  and  UACC62)  (a,  c,  e),  and  the  murine  melanoma  cell  line 
YUMM  1.7  (g)  and  corresponding  vemurafenib-resistant  derivatives  (A375R, 
Colo800R,  UACC62r  and  YUMM1.7R)  under  increasing  concentrations  of 
vemurafenib.  Immunoblotting  of  phosphorylated  ERK  protein  levels  in 
indicated  melanoma  cell  lines  in  the  presence  of  increasing  concentrations  of 
vemurafenib  (b,  d,  f,  h).  i,  Immunoblotting  of  phosphorylated  ERK  and 


phosphorylated  AKTS473  protein  levels  in  HCC827  lung  adenocarcinoma  cells 
in  the  presence  of  increasing  concentrations  of  erlotinib.  j,  Immunoblotting  of 
phosphorylated  ERK  protein  levels  in  H3122  lung  adenocarcinoma  cells  in 
the  presence  of  increasing  concentrations  of  crizotinib.  k,  Immunoblotting  of 
protein  levels  of  MET,  EGFR,  BRAF,  PDGFRb,  phosphorylated  AKT  and 
phosphorylated  ERK  in  vemurafenib -sensitive  and  -resistant  pairs  of  human 
melanoma  cell  lines  (A375,  C0I08OO  and  UACC62). 
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Extended  Data  Table  1  |  Clinical  data  for  tissue  donor  subjects 


Study  site 

Pt# 

Bx  samples 

Age  &  Sex 

Stage 

Dose  (mg) 

BOR 

PFS  (days) 

Bx  site 

UCLA 

1  TG 

Baseline 

Day  7 

51 M 

Mic 

960  bid  vemurafenib 

-21% 

108 

SC,  scalp 

SC,  scalp 

2  JCC 

Baseline 

Day  15 

44M 

Mic 

960  bid  vemurafenib 
+60  qd  cobimetinib 

-63% 

Current  response 

SC,  abdomen 

SC,  abdomen 

3  YAU 

Baseline 

Day  22 

26F 

Mic 

960  bid  vemurafenib 
+60  qd  cobimetinib 

-46% 

145 

Dermal/SC,  abdomen 
Dermal/SC,  clavicle 

bid,  twice  daily;  BOR,  best  overall  response;  Bx,  biopsy;  F,  female;  M,  male;  PFS,  progression-free  survival;  Pt,  patient;  qd,  daily;  SC,  subcutaneous;  UCLA,  University  of  California,  Los  Angeles. 


©2015  Macmillan  Publishers  Limited.  All  rights  reserved 


